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Abstract 


This report describes the results of a research study of the feasibility 
of conducting experiments in the Spacelab on the ignition of, and flame 
spread over, liquid fuel pools which are initially at a temperature 
lower than the fuel's flash point temperature. Theories are developed 
for the ignition and flame spread processes and experiments are conducted 
to understand the factors influencing the ignition process and the spread 
rate. The results are employed to devise a conceptual Spacelab experiment 
which is expected to be feasible for a safe conduct and to be suitable for 
obtaining crucial data on the concerned processes. 
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I. Introduction 


The research reported here deals with the issues of feasibility of 
studying the burning of liquid fuel pools in reduced gravity conditions 
facilitated by the Spacelab. This research was sponsored by NASA under 
contract No. NAS3-21018 to the University of Notre Dame, Department of 
Aerospace and Mechanical Engineering. Based on a previous study [1] 
(under NASA Contract No. NAS3-20G87, NASA Rept. No. CR-135234, Notre Dame 
TR. No. 77-33, June 1977) the specific topics of current study evolved to 
be ignition of, and flame spread over, fuel pools whose temperature is 
below the fuel’s flash point. 


1. Ignition : If a heated rod or a well-defined flame, or any other heat 

source were placed parallel to, and in the proximity of, the surface of a 
fuel at a temperature lower than its flash point temperature, a transient 
free convective flow is set-up around this source in the gas phase. This 
flow is strongly dependent upon gravity and establishes a thermal communi- 
cation between the source and the fuel surface if the source is suffi- 
ciently close to the surface. The heat flux to and temperature of the 
fuel surface will be maximum directly under the source, and they will 
decay monotonically as one marches away from the source in a direction 
normal to its axis. Due to the temperature gradients along the liquid 
surface, surface tension gradients are set-up to produce a flow in the 
condensed phase. Buoyancy, shear, and inertia also come into play in 
this flow which forms two symmetric cylindrical cells as shown in Figure 
1 . 

Given this situation, to determine the flow patterns, heat transfer 
distribution, resultant fuel vaporization, vapor dispersion in the gas 
phase, and ignition of the mixture constitute problems of extreme im- 
portance both in science and practice. In the present project the 
development of flow and heat transfer is studied to examine the influence 
of the properties of the fuel, fuel layer thickness, gravity, heat source 
strength and dimensions, etc. on the time to ignition. Experiments are 
discussed in Chapter V and Appendix E whereas the results of a theory 
presented in Appendix D are summarized in Chapter IV. 

2. Flame Propagation : Once a liquid fuel is ignited as above, say at 

one end of a long tray, the flame would propagate over the liquid surface. 

Flame propagation over a horizontal pool of fuel follows one of two 
possible distinct patterns [5] depending upon whether the fuel’s initial 
temperature is above or below its flash point. When the temperature is 
above the flash point a combustible gas mixture exists above the liquid 
surface. After ignition a flame would then propagate through the 
combustible mixture under the control of the gas phase mechanisms of 
heat and mass transfer. 
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When the fuel temperature is below the flash point, however, the 
propagating flame will have to preheat the fuel ahead of it continuously, 
this preheating being just adequate to raise the adjacently located fuel 
element to the flash point. The flame then creeps forward to this heated 
element and proceeds to heat the next element. Keeping aside the impor- 
tance of thermal radiation, this preheating can occur by conduction and 
convection in both gas and liquid phases. The gas-phase heat transfer 
becomes inconsequential since the gas-phase conductivity is relatively 
negligible and since the gas-phase convective flows oppose the flame 
propagation direction. The liquid-phase conduction is smaller by about 
an order of magnitude [2] than the liquid-phase convection which occurs 
in the seme direction as flame propagation 

Experiments of Roberts [10,11], and Mackinven [6] showed that the 
condensed phase flow is induced by the flame itself and that the two of 
the most influential parameters are the liquid viscosity and the sensi- 
tivity of liquid surface-tension to the surface temperature. The liquid 
flow patterns set up by the combined influence of gravitation and surface 
tension assume the primary responsibility to transfer energy ahead of the 
flame front. A qualitative sketch of this flow is presented in Figure 2. 
Experiments on flame spread are discussed in Chapter V and Appendix E 
whereas the essential results of a theory developed in Appendix C are 
presented in Chapter III* 

3. Justification : The Spacelab experimentation on liquid fuel ignition 

and flame propagation promises to yield valuable, otherwise unattainable, 
data. As discussed below, the heat and mass transfer processes leading 
to ignition and flame spread involve several competing or cooperative 
forces. Inasmuch as the total problem is much too complex, simplifications 
are required to be made in order to develop theories of ignition and of 
flame spread. These simplifications are only possible upon the knowledge 
of the relative importance of the various forces. This knowledge is re- 
quired also to test the validity of assumptions made in theoretical 
studies. Earth-based experiments are unsuitable or inconvenient in 
attempts to isolate these forces and understand their role. Thus it is 
expectable that the Spacelab data are crucial in advancing the state-of- 
the-art of the liquid fuel combustion problem. 

a. Theoretical : The unusual motions of the liquid resulting from 

an intricate interaction of surface-tension forces with gravitational, 
viscous, and inertial forces are of interest not only in science but also 
in practical applications. Conduct of liquid flow experiments iu space 
gives an extra degree of variability of the gravitational constant to 
control the Froude number (i.e., inertial forces/gravitational forces). 

Bond number (i.e., gravitational forces/surf ace-tension forces), and no- 
name number (i.e., gravitational forces/viscous forces) [1]. The data 
obtained in the Spacelab over extended ranges of these dimensionless 
quantities will offer a powerful basis to verify various theoretical 
hypotheses on the special liquid motions. 

Over and beyond aiding to satisfy one’s scientific curiosity about 
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these motions, the Spacelab experiments are expected to lead also to a 
better understanding of a host of practical situations involving liquids, 
slurries and emulsions subjected to heating or cooling, and evaporating 
with or without chemical reactions. These chemical reactions may pertain 
to combustion or to any of the numerous chemical engineering processes. 

Both ignition of liquid fuel pools and flame spread over them are 
expected to involve important transient effects of heat, mass and momentum 
transfer. These transient effects are strongly dependent on the 
gravitational level. At high values of g, it is also expectable that 
three-dimensional phenomena would overwhelm. Conduct of experiment in 
space where reduced gravity is possible is anticipated to enable one to 
delineate these and other complicating secondary aspects of the studied 
processes of ignition and flame spread over liquid surfaces. 

Classman's work [5-7] pointed out that surface-tension-driven flows 
in the liquid are of paramount importance in the convective heat transfer 
responsible for the flame propagation. Once again, the liquid viscosity, 
fuel layer depth, sensitivity of surface tension to the surface temperature 
are some of the critical factors in this heat transfer process. 

Besides Glassman's group, Torrance [8,9] and Roberts [10,11] are among 
the authors who studied surface-tension flows. The available analyses 
are, however, unsuitable to predict the role of gravity in either aiding 
or opposing the surface-tension-driven flow. Conduct of experiments 
under low gravity conditions available in the Spacelab will permit the 
use of conveniently deeper pools to study the shallowness effects, for 
then the pools are deep enough to be instrumented or visualized but at 
the same time the low gravity preserves the physics by keeping low the 
governing Grashof, Reynolds and Marangoni numbers. (These nondimensional 
numbers are defined in Table I). 

Torrance [8,9] computed the flame spread Reynolds number as a 
function of Marangoni and Grashof numbers. His main finding was that 
buoyancy and surface-tension work hand-in-hand to produce a thermal con- 
vection cell which enables the flame to propagate forward. The effect 
of buoyancy is small and nearly indistinguishable for Gr * 10 . At 
Grashof number of about 10^, buoyancy appears to influence the flows 
about as much as surface-tension. In contrast to this, Glassman [7] 
assumes gravity's influence is realized through hydrostatic pressure 
which opposes the surface-tension-driven flow. Under these circumstances 
the absence of gravity would significantly affect the flame spread rate. 

It is clear that there is not a firm understanding of the role gravity 
plays in flame spreading processes. Carefully conducted experiments in 
reduced gravity offers the prospect of shedding light on these complex 
phenomena , 

Our considerable search into the chemistry, physics and engineering 
literature pointed out only a few references to the ignition problem 
[2-4] . Dimensional and hueristic analyses point out the role played by 
gravity [1], but to quantify the expected results, experiments in the 
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prolonged programmed gravity environment of the Spacelab are highly suit- 
able. 

Historically, the science of combustion possessed a special character 
in that the incubation period is relatively short for scientific concepts 
to become translated into important practical applications. Thus, progress 
made in combustion science via Spacelab experimentation may be expected 
to lead to prominent technological developments. 

a. Experimental : If variable or low gravity experiments are thus 
deemed justifyable on the basis of theoretical arguments, one might 
question why one can not simulate the reduced gravity conditions in earth- 
based laboratories; and why one would require experimentation in the 
Spacelab. In Appendix A, a variety of currently used techniques of simu- 
lating low gravity conditions are reviewed and their various features are 
discussed. Some of these techniques are indeed suitable for studying 
some short-lived processes of fluid physics and heat transfer. However, 
if the process under study involves characteristic time scales of the 
order of minutes as do the liquid fuel ignition and flame spread, the 
state-of-the art simulation techniques become unsuitable, for the low-g 
duration offered by them amounts to the order of 10 seconds or so. Thus 
it appears to be eminently justified on the basis of experimental suit- 
ability to advocate liquid fuel pool experiments in the Spacelab. 

b. General : Over and above these preceding analytical and experi- 

mental justifications of the study of fuel liquid ignition and flame 
spread in the Spacelab, there is also a more general justification. 

This pertains to the need to establish fire safety in the low-gravity 
conditions of our space vehicles in general, and in the Space Shuttle in 
particular. These various installations are essential and extremely 
expensive. Their cargo is ultra-sophisticated scientific equipment 
whose material value is greatly overshadowed by the value of the critical 
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data and information they bear to generate towards the advancement of 
science and solution of technological problems f of mankind. Even more 
important is the life safety of the Spacelab personnel. Intrinsically, 
the options available to fight an accidental fire, or to egress to safety 
from it, are limited. Thus, again, it is critical to gain a thorough 
understanding of the processes of fire initiation and growth at low gravity 
conditions. 

Thus it is clear that the Spacelab experiments on liquid fuel ignition 
and flame propagation are desirable. The theories developed in Appendices 
C and D and discussed in Chapter III and IV as well as the experiments 
described in Chapter V and Appendix E point out some further reasons why 
the prolonged conditions of low- or zero-gravity are desired for the study 
of ignition and flame spread. 



II. Feasibility Issues 


Once the need for further study of the liquid fuel ignition and flame 
spread processes is substantiated [1] and the conduct of the experiment in 
the essential environment provided by the Spacelab is justified analyti- 
cally and experimentally, our immediate next concern is to establish the 
feasibility of the Spacelab experiment, 

Typical questions related to feasibility examination are as follows: 
which particular fuels will yield the most useful results? Under what 
particular ranges of the independently variable conditions do these useful 
results arise? How safe is the conduct of the proposed experiment in the 
Spacelab? How large will the fuel pool have to be to give meaningful 
results? How would the fuel be transported, stored and loaded into the 
pool tray? How full would the tray be filled; is a 100% full tray possible 
under the influence of the curvature of surface resulting from surface- 
tension? How will the flow patterns be observed and recorded in both the 
liquid and gas phases? How are the temperature, velocity and concentration 
fields to be quantitatively measured in the Spacelab experiment? Are 
there any physical features of the ignition and flame spread processes 
which can be studied by observing the thermo-fluid -dynamic response of 
inert liquids, say water? What exactly are the hypothetical mechanisms 
of ignition and flame spread of liquid fuels? 

The paramount issue in dealing with liquid fuels in the Spacelab 
appears to be that concerned with safety. It infringes upon decisions 
related to nearly all other design aspects of the experiments. 

1. Choice of Fuels : The factors involved in this choice are: the 

flash point of the fuel in air; the fuel's corrosivity and toxicity; its 
known thermophysical properties; etc. The table of properties summarized 
in Appendix B leads one to believe that such fuels as benzene, xylene and 
toluene are undesirable because of their extreme reactivity and 
corrosiveness. Our experiments, described elsewhere in this report, show 
that alcohols are suitable, especially butanol. 

2. Pool Dimensions : The fuel pan dimensions have to be small enough 

to alleviate a catastrophe in the event of an accidental spill and yet 
large enough to minimize two-dimensional effects of the side walls. 
Glassman’s work [6] indicates that not unless a width of about 0.60 meters 
or more is employed will one realize a flame spread process uninfluenced 
by the side-walls. Since this magnitude of a width is clearly unaccept- 
able for the Spacelab experiment, we propose to use a small width but by 
developing theoretically and verifying experimentally at 1-g the influence 
of pool width. A three-dimensional theory has to be developed to 
segregate clearly the effects of the pool width on the ignition and 
flame spread processes. Our feasibility study experiments on butanol 
pools ranged in width from about 2 cm to 10 cm, as reported in Appendix 

E of this report. 
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The depth of the pool has to be so chosen as to avoid impeding the 
surface-tension-driven flow cell whose dimension depends on various 
properties of the fuel and flame as discussed in Appendix C. For butanol 
this depth was found to be about 1 to 2 cm from our experiments. 

Finally, the length of the tray has to permit observations free of 
the region in which the ignition transients dominate and from the region 
where the upstream boundary effects affect the condensed phase flow field. 
An overall length of 30-35 cm appears to satisfy these restraints for 
butanol pools. 


3. Fuel Transportation, Storage and Handling : Based on the Claussius 
Clapeoron relation between the saturation vapor pressure and temperature, 
it appears that the fuel could be transported and stored in a variable 
pressure container. The flash point temperature is merely the temperature 
of the liquid at which, corresponding to the prescribed total pressure, a 
critical vapor pressure is developed over the liquid surface. The 
C.riticality is stipulated by the lean limit of flammability. It appears 
that for a given transportation and storage temperature of the liquid, 

the lower the maintained total pressure the higher will be the flash point 
relative to the storage temperature. Thus one can transport and store 
fuel at a low pressure safely and then introduce it into the pan at a 
higher pressure, the level of which is determined by the desired flash 
point relative to the supply temperature. This technique of transporting 
the fuel safely is suitable for fuels whose normal flash point is close 
to the ambient temperature of the experimental facility. For fuels such 
as butanol whose flash point is considerably higher than the supply or 
ambient temperature, safe transport and storage are much easier. 

A method of handling the fuel in this safer situation is described 
in the conceptual design of Spacelab experiment. This method recognizes 
that the volume of liquid required per experiment in the chosen pan is 
about 270 m£. This volume corresponds to a sphere of diameter * 8 cm. If 
a pair of hemispherical clamshells are fabricated with flanges which can 
be bolted together to form a spherical chamber with an elastic membrane 
in-between them, the required amount of fuel could be loaded on one side 
of this membrane. Employing compressed air, the fuel could be displaced 
into the pool tray or sucked back into the sphere after an experiment. 

The technique definitely offers a controlled and safe transportation, 
storage, and handling of fresh fuel and stowage of used fuel. 

4. The Filling Process, Stability of the Free Surface : With the fuel 
transported/stored employing the spherical container discussed above, the 
filling process is quite simple. As the pan is filled, What are the 
characteristics of the free surface as to its stability? In attempts 

to answer this question, Labus at NASA [12] has conducted a dimensional 
analysis and some simple experiments. Liquid was introduced into a 
precisely machined stainless steel pan which was then tipped slowly 
until the fluid flowed out. Denoting, for convenience, the plane of the 



undisturbed surface to be 'horizontal' and its normal to be 'vertical', 
the limiting tilt angles 0 are correlated on a horizontal Bond number 
(g sin0 p/cr) versus vertical Bond number (g cos 0.R 2 p/a) plane 
(where g is gravitation or gravitational jitter, 0 is angle of tilting 
from horizontal, R is pan radius, p is liquid density and a is surface- 
tension), The results of the experiments showed that in the limit of a 
very small vertical acceleration environment, the critical horizontal 
Bond number required to retain the fluid in the pan has to be less than 
2 to 3. For butanol the density p s 0.81 g/cm^ and the surface-tension 
a * 25 dynes/cm. For a pool of width 2.54 cm, the half -width R = 1.27 cm. 
Thus if the critical horizontal Bond number of 2 to 3 observed by Labus 
is valid for butanol , the maximum permitted gravity and inclination to 
avoid spilling by width-wise 'rocking' of the pool are given by gsin0 * 2 
to 3 2 times the quantity (a/pR 2 ); (8 sin 9) raax>perin;Ltted - 38.27 - 57.4 
cm/s . Noting that smaller values of either gravity or g-jitter and of 
inclination would ensure containment of the liquid in the pan, a factor 
of safety of 2 would result in a maximum permitted (gsin0) of about 19.13 
to 28.7 cm/s . A horizontal pan thus could be subjected to a g-jitter 
of as much as 0.02 to 0.03 times the normal gravity of 981 cm/s 2 before 
causing a spill. The g-jitter values expected in the Spacelab due to 
extraneous disturbances are of the order 10“^ times the normal g. The 
pan width of 2.54 cm thus appears to provide safety from spillage due to 
g-jitter. 

If the same arguments are made based on the length-wise 'rocking' 
of the pan of length 2R=35 cm, the maximum permitted value of (gsin0) 
to avoid spillage is about 0.10 to 0.15 cm/s^. The maximum permissible 
g-jitter to avoid spillage due to length-wise disturbances is thus about 
1 x 10~4 to 1.5x10"^ times the normal gravity of 981 cm/s 2 . Thus the 
spillage is a couple of orders of magnitude easier to occur due to length- 
wise disturbances than due to width-wise disturbances.. Even the more 
vulnerable length-wise spillage will not be a problem with a 35 cm long 
butanol pool if the g-jitter is 10~5g o Q r lower. 

5. Methods of Flow Visualization and Measurement : The current state- 

of-the-art of flow, and heat and mass transfer measurement techniques are 
reviewed by this author in reference 13. In the Spacelab experiment, we 
expect to employ particle trace techniques of flow visualization. The 
tracer, most probably metallic particles or aluminum oxide flakes, will 
be premixed in the fuel container. 
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III. Flame Propagation Over Liquid Pools 


The basic physical picture of the flame propagation over a liquid 
fuel layer whose temperature is below its flash point is shown in Figure 
2. The pool is assumed to be infinite in extent normal to the plane of 
the paper, i.e., the pool width is assumed to be irrelevant. The depth 
of the undisturbed pool is h^ so that y*0 forms the rigid base of the 
pool and y*h fo refers to the liquid-air interface. Suppose a flame is 
initiated at x = - °° and is permitted to steadily and uniformly propagate 
in the +x direction at a speed of U. If the frame of reference x=0 is 
fixed on the spreading flame front, one encounters a steady phenomenon 
in which the flame front is at rest while the fuel and gas-phase are fed 
backwards into the flame. 

The surface temperature of the fuel layer decreases rapidly with 
distance upstream of the flame front; at the flame front x®0 this 
temperature will be equal to the fuel's flash point whereas at some 
positive x-value the temperature will be equal to the lower supply 
temperature. Since the surface-tension of most liquid fuels decreases 
with increasing temperatures, the nonunifomity in surface temperature 
would result in a nonuniformity in surface-tension. Thus the effect of 
surface-tension gradient in the x-direction is essentially to pull the 
surface liquid away from the flame front. As the hot liquid is thus 
carried forward in the direction of flame propagation, energy is con- 
vected to continue the propagation further. 

Sirignano and Glassman [7] and Torrance [8,9] carried out theoretical 
an 3ysis of the flame propagation controlled by surface-tension flow. 

There are four forces of interest in the problem; these are: inertial, 

gravitational, viscous, and surface-tensional. These four forces result 
in three independent parameters to describe the problem. For flows driven 
by surface-tension, the order of magnitude of the velocity is determined 
by a balance between the surface-tension force and the viscous force at 
the surface. Thus the ratio of: surface-tension force to viscous force is 
of the order of unity. As shown in Table I, the ratio of inertia force 
to viscous force is known as Reynolds number; the ratio of buoyancy to 
inertial force is Froude number; the product of the ratio of buoyancy to 
viscous force and the ratio of inertia to viscous force is known as 
Grashof number; and finally the ratio of heat transfer due to surface- 
tension-driven convection to conductive heat flux is known as Marangoni 
number. 

Torrance [8,9] compared the flame spread Reynolds number as a function 
of Marangoni and Grashof numbers. His main finding was that buoyancy and 
surface-tension work hand-in-hand to produce a thermal convection cell 
which enables the flame to propagate forward. The effect of buoyancy is 
small and nearly indistinguishable for Gr - 10^ . At a Grashof number of 
about 10^, buoyancy appears to influence the flows about as much as sur- 
face tension. He also predicted that the surface temperature would 
decrease with distance ahead of the flame front in an almost linear 
fashion and would become uniform at larger distances. The characteristic 


PACE tP INT^YIGHALLY Sli® 


11 



length for this surface temperature, is independent of the liquid layer 
depth but is inversely proportional to the flame spread speed. In 
contrast to this, Sirignano [7] assumes that gravity's influence is 
realized through hydrostatic pressure which opposes the surface-tension- 
driven flow. Under these circumstances the absence of gravity would 
significantly affect the flame spread rate. 

In Appendix C a theoretical model of the surface-tension-controlled 
flame spread over a liauid fuel surface is presented. Sirignano [7] 
developed a similar theory but by ignoring buoyancy and by focusing on 
the hydrodynamic part of the problem. We considered, following Birikh 
[47], the energy conservation equation to obtain the temperature field 
from which we can account for the buoyancy term in the momentum equation. 
Torrance [8,9] studied this problem comprehensively by employing numeri- 
cal methods of solution. Even though our solution pertains only to the 
field ahead of the steadily propagating flame, it is of a closed-form. 

The essence of Appendix C for the present feasibility study is repre- 
sented by the following results of nondimens ional velocity and temperature 
fields. 


;=-Re+Ci[^ + 2 

„ 2 

- PrG l T Cr - ~ 1 

T= (1 + C 1 x)y+ a(y) + * b(y)J 


a(y) 5 (100 y 7 - 567 y 5 + 490 y 4 + 99 y 3 - 122 y) 
b(y) = (9 y 5 - 10 y 4 - 3 y 3 + 4 y) 

where the nondimens ional quantities listed in Table II are defined as 
below: 

u = ph u/y nondimensional velocity 

flame spread 
temperature 
Grashof number 
surface-tension parameter 


Pr = „c p / K 


Re 5 ph^ u/y 
T 5 (T-TJ/(T 0 -TJ 
Gr e gh 3 p 2 6(T -T )/y 2 

00 O co 

* = * (I 0 -TJ ph„/w 2 
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= x/h^ (0 1. x <_ - 1/ C^) 

y = y/hoo (o < y < 1) 

C 1 5 « Re/3 (i --$g Gr) 

u is the x-directional velocity, U is the flame spread speed* h^ is the 
undisturbed fuel layer thickness, p is fuel density; p is its dynamic 

viscosity; I 0 is its flash point; 3 is its volumetric expansion 

coefficient; <f> is the sensitivity of surface-tension a to temperature, 

<J> = dcr/dT, a negative quantity; Pr is the liquid's Prandtl number; K. is 
thermal conductivity; and Cp is specific heat, g is gravitational 
constant and To, is fuel supply temperature. 

The nondimensional surface-tension parameter <j> is related to the 
Marangoni number mentioned earlier in a straight-forward way. To obtain 
the constant Cp in the form given above constitutes the essence of the 
theory presented in Appendix C. For x > l/C-^, u •- -Re and T = 0 uniformly. 

Figure 3 shows the velocity profile in the fuel layer as predicted. 
Inasmuch as the temperature profile is linear with respect to x, this 
velocity profile is independent of x in the domain 0 £ x <_ - 1/C^. The 

flow reversal is evident from Fig. 3. Additionally it can be seen that 

the effects of increased surface tension sensitivity and reduced gravity 
are cooperative. The parameter Gr/ij> in this figure is a measure of the 
relative intensity of gravitational and surface-tension forces. 

Gr/<(> = 0 refers to zero-gravity situation while Gr/ <j> = - °° refers to 
zero-surface-tension-gradient situation. Gr/cj> = -60 approximately 
corresponds to the situation encountered in the earth-based laboratory 
With a 2.5 cm thick layer of a fuel such as butanol. 

Reduction or absence of gravity, in general increases the flow velo- 
city at the surface. This increase, according to the continuity equation, 
is only possible at the expense of a decrease in maximum negative velocity 
at a y of about 1/3. Note also that the velocity gradient 9u/3y at the 
surface approaches zero as gravity increases towards infinity relative to 
the surf ace- tension-grad lent. This is understandable, for any shear at 
the free surface has to be balanced by the surface-tension force. It is 
interesting to note that our resultant velocity profile, when Gr/lj> is set 
to zero, coincides with the Levich profile. And finally, the influence of 
reduced gravity on the velocity profile is rather abrupt; this can be 
seen by noticing that whereas the difference between the velocity profiles 
of Gr/$ = 0 and Gr/cp = -60 is considerable, the difference between the 
profiles of Gr/(j> = -60 and - °° is relatively minor. 

The temperature profiles are shown in Figure 4 (flame spread Reynolds 
number is 2) and 5 (Re =10). The surface-tension parameter $ and Prandtl 
number Pr are kept fixed respectively at -100 and 45 but three different 
Grashof numbers are considered. The temperature, of course, varies 
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FIGURE 3: THE VELOCITY PROFILE IN THE FUEL LAYER AS INFLUENCED BY GRAVITY, SURFACE TENSION 

AND FLAME SPEED. 



FIGURE 4: TEMPERATURE PROFILES IN THE FUEL LAYER AS INFLUENCED BY 

GRAVITY AT FIXED Re, $ and Pr. 
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between a maximum of unity and a minimum of zero. 

It is clear that the larger the gravity the greater the penetration 
of the thermal disturbance into the liquid layer. It is also clear that 
slowly propagating flames will result in a deeper penetration of the 
thermal effects. Furthermore, the temperature effects reach deeper as 
one approaches (in x-direction) the flame’s leading edge. 

In Figure 6 the surface temperature as a function of x, distance from 
the flame front, is shown specifically as dependent upon the flame spread 
Reynolds number, surface-tension parameter and Grashof number. The 
linear relation is an assumption in the present theory. However simpler 
analyses lead to the expectation of this linearity even though one’s 
intuition can not readily lead to it. The linearity is altogether con- 
sistent with Torrance’s [8,9] numerical solutions of the full governing 
equations. 

To find the slope of the T s (x) linear relation constitutes the 
essence of the theory presented in Appendix C. This slope, denoted by 
Cj_, is given above. Faster propagating flames, over pools of fuels with 
lower surface-tension parameter $, at lower levels of gravity produce 
steeper decay of the Surface-temperature with distance x. The influence 
of the gravity and flame spread rate are much stronger on this steepness 
than the influence of surface-tension. 

The slope has the importance that it indicates the. extent of for- 
ward thermal influence, x oo =-l/C^, exerted by the combined Gr and cf> 
effects. It is thus useful in delineating the length required of the 
Spacelab pool experiment. For alcohols, $ =-100 to -200 and Re 8 1 to 
10. Earth-based experiment on pools, about 2 to 3 cm deep, of alcohol 
have a Gr of about 6000. Thus one can see from Figure 6 that the region 
of forward thermal influence is about an order of magnitude shorter in 
zero-g than in normal-g. 

The x to - dependency on Grashof and Reynolds numbers, and the surface- 
tension parameter is portrayed in Figures 7 and 8 in order to graphically 
show the involved sensitivities. We must stress again that the influence 
of reduced gravity appears to be to reduce the length of thermally dis- 
turbed region. Reduction in surface-tension gradients and increase in 
flame spread velocity also have the same consequence albeit less strong. 

How do we employ these inferences in choosing the dimensions of the 
Spacelab experiment? The viscous forees at the base of the pool tend to 
elongate the recirculatory mass of the fluid. For any chosen depth, the 
distance ahead of the flame front to which the motion and thermal distur- 
bances extend is of the order of ten times the depth in an earth-based 
laboratory whereas it is of the order of one to three times the depth at 
zero-gravity. Thus choosing a pool of depth 2.54 cm (arbitrarily picked 
for convenience) the region of length-wise interest is expected to be 
about 2.54 to 7.5 cm at zero-gravity. For suitable observations, there- 
fore, a 15-20 cm observation length of the pool is deemed adequate with 
about 7-15 cm ignition transient length and about 7.5 cm length for 
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alleviating the suppression effects of the upstream boundary. Ignition 
transient length is the distance from the ignition point in which the 
flame would induce itself to a state of steady propagation. From the 
analysis of the ignition problem presented in Appendix D and summarized 
in Chapter IV and from the experiments discussed in Chapter V and in 
Appendix E, this length is expected to be about 7-15 cm in zero-gravity 
and about 10-15 cm at normal gravity* 

The Spacelab pool, according to these considerations, is to be about 
30-40 cm long, 2.54 cm deep and an arbitrary width of 2.54 cm. 



IV. Ignition 


As a quiescent fuel layer of thickness h^ is subjected to a heat 
source of strength q' watts per unit length in the pool-width direction, 
heat transfer occurs to the liquid transiently, surface-tension-driven 
flow would ensue along with natural convection and there would arise an 
instant at which the liquid is locally raised to its flash point to re- 
sult in ignition. A mathematical model is presented for this problem in 
Appendix D. The main results of this theory pertain to the transient 
development of the extent of flow and thermal perturbation, the surface 
velocity and the liquid temperature averaged over the region of distur- 
bance. The following results are obtained. 

Disturbance size; 



Surface velocity: 



Liquid temperature at surface: 



d is the diameter of ignition source, <j> is surface-tension sensitivity 
to temperature, Gr is Grashof number and q is power fed to the ignitor, 
all dimensionless as def ined_below along with the nondimensional time 
t, dimension of disturbance A and velocity u. 

x 5 x/h^ ; d = d/h TO ; A s- A/h m 

u = ph u J]i ; t 5 pt/ph ; 

s 00 s 00 

+ s ♦Cr o -T„> Ph> 2 S ; = «47» C p (T o -TJ; 

Gr = gh 3 p 2 S (I -T )/u 2 . 

cp o o° 

x is distance along the liquid surface; d* is essentially the diameter 
of the ignition source; h^ is the pool depth, A is cell dimension; 



Ug is surface velocity; p is liquid density; y is dynamic viscosity; t 

is time; <f> is da/dT, a being surface-tension; T Q is fuel flash point 
temperature; T ro is fuel supply temperature; q' is power input to the 
ignitor per unit length in the direction of the pool's width; f is 
fraction of q' delivered to the pool; C p is heat capacity of the liquid; 
g is gravitational acceleration, and 0 is the volumetric expansion 
coefficient of the liquid. It is interesting that the dimension of the 
disturbance would grow with time raised to the power 2/5. The greater 
the Grashof number and - <£ , the larger the growth. A heat source large 
in size and strength would also enhance the growth. If <J> = -100, the A 
at the same heating rate and time at zero-g will be about 75% of its 
value at the earth's gravity (i.e. , Gr* 6000). If $ = - 200, this fraction 
becomes about 83%. This somewhat unexpectedly weak influence of gravity 
mnkes the earth-based studies of preignition flow patterns useful in 
interpreting the results of the Spacelab experimental data. 

Note from the surface velocity equation given above that at time near 
zero, A being small, an impulsive motion is set up in the liquid surface. 
This prediction is consistent with the findings of MacDonald [4]. 

As the transient build-up of motion and heating continues, ignition 
of the liquid is said to have occurred at the instant when the surface 
temperature at x = 0 reaches the flash point. Thus setting T g = 1 at 
x = 0 in the solution, the ignition conditions (denoted by asterisk) are 
obtained . 



For alcohols - 100 to -200. If q’ lies between 0.1 and 1.0 kW/m, 
q z 1.5 - 15 with alcohol supplied at about 300 °K. The heat wire dia- 
meter can be ranged such that d lies between 0.05 and 0.20. The gravity 
and surface-tension effects conspire together to determine the ignition 
time and the extent of the disturbed region as shown in Figure 9. 

The ignition time and extent of disturbed fluid are larger with a 
thicker ignitor, keeping all else fixed. The surface velocity at this 
time directly under the heater, however, will be lower with thicker 
ignitor. Furthermore, the ignition time is inversely proportional to the 
cube of heating rate. This strong dependency arises because the fluid 
velocity itself is a linearly increasing function of q; the resultant 
enhancement in the convective dissipation of energy is the reason why it 



takes much longer to arrive at the ignition condition. 

Note from Fig. 9 that ignition time as well as the flow region 
dimension are significantly reduced by reducing gravity to zero. For 
example, with a of 15 and d of 0.20, t* *0.4 in an earth-based 
laboratory and 0.025 in the Spacelab; the corresponding values of A* are 
respectively 2.6 and 0.67. At other values of q and d similar differences 
exist between the earth-based laboratory and the Spacelab as listed below: 


d 

d 

earth 

t* 

Spacelab 

earth 

A* 

Spacelab 

1.5 

0.04 

16 

1 

5.4 


1.3 

1.5 

0,20 

400 

25 

26 


6.5 

15 

0.04 

0.016 

0.001 

0.54 


0.13 

15 

0.20 

0.4 

0.025 

2.6 


0.65 


From this table, one can see that under expectable heater and fuel layer 
circumstances, a maximum. A* of about 6 is realized in a reduced gravity 
environment. Thus, if our Spacelab pool is arbitrarily chosen to be 
2.54 cm deep, for convenience, the preignition disturbance will extend 
as far as 0.33 cm to 15 cm. The ignition time t* is expected to lie 
in the range 0.125 sec. to several minutes. 
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V. Experimental 


The experimental part of this project, described in detail in 
Appendix E, led to the following conclusions of relevance in the 
feasibility issues. 

(a) A tray of dimensions 2.54 cm wide, 2.54 cm deep and about 35 cm 
long appears to be suitable to obtain meaningful data on preignition and 
flame spread. The depth is chosen arbitrarily for convenience in handling 
the fuel. The width effects appear to be amenable for isolation by a 

3-d theoretical study and earth-based experiments. 

(b) Disturbances in the fluid during the preignition heating appear 
to extend to about 10 cm. 

(c) The steadily propagating flame on 2.54 cm deep butanol appears 
to produce flow and thermal perturbations to a distance of about 10-15 
cm ahead of the flame front. This distance is in reasonable agreement 
with the flame spread theory. In zero-gravity environment, the theory 
leads us to expect that this distance will be much smaller, i.e., about 
1.5 - 7 .5 cm. 

(d) The tray could be fabricated out of an aluminum frame with in- 
serted plateglass sidewalls. 

(e) The environmental temperature and humidity appear to be of little 
consequence on the observed flame spread rate. The room drafts, however, 
are expected to disturb the propagating flame. 

(f) The nature of fuel effects the flame spread -related and ignition- 
related processes drastically. Whereas the flash point temperature at the 
operational air pressure is the most important property of the fuel, the 
liquid-air interfacial tension and its sensitivity to temperature are 
also crucial. Other properties of concern include the liquid density, 
specific heat, thermal conductivity, viscosity, etc. 

(g) Xylene, a ring compound is found in our experiments to be 
extremely corrosive. 

(h) Commercial kerosene, a usually poorly defined mixture of several 
hydrocarbons, is quite suitable for preignition and flame spread experi- 
ments, but the data can not be employed to test theoretical hypothesis 
because of the unknown physico-chemical properties. 

(i) Butanol is found to be a well-defined and suitable fuel for 
further work. The only minor drawback of this fuel is its ability to 
absorb and dissolve water vapor from the flames. 

(j) The flame spread rates are erratically large as the initial 
supply temperature of the fuel approaches its flash point. The most 
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orderly observation of the butanol flame spread process in 1 atm air can 
be done at temperatures below about 38°C. 

(k) The free-board height is to be chosen either to be zero (i.e., 
the fuel surface is flush with the top edge of the tray walls) or to be 
larger than 6 mm (i.e., the fuel surface is over 6 mm below the top edge 
of the tray walls). 

(l) Flame spread rate is larger in trays of larger widths. The side 
wail effects appear to be greater at larger fuel supply temperatures. 

A tray width of 2.54 cm appears to be suitable for the Spacelab experi- 
ments. The observed flame spread rates and flow patterns could then be 
corrected for the width-effects. 

(m) There appears to be a linear relationship between the flame spread 
rate and the pool depth. The slope is larger for larger fuel supply 
temperature. For pool depths grenter than about 1.5 cm the linearity is 
well-established. 

(n) The flame spread rate may be influenced by the ignitor disturbances 
to an extent of about 10-15 cm on the ignitor end. However, the effect of 
the upstream end-wall (on the opposite end) may be felt by the flame for 
about 10-15 cm also. A length of 15 cm in the middle of a 35 cm long tray 
is deemed long enough to obtain data on a flame whose spread rate is 
uninfluenced by the end effects. 

(o) Butanol fuel exhibits quite vivid liquid flow patterns when ex- 
posed to a small propane flame (about 5 mm in diameter and 1.5 cm in 
length) placed at a distance of about 5 mm above the surface. Observation 
of the transient development of the surface-tension-driven convection 
flow is possible and the ignition time is of the order of 0.1 to 1 sec. 

(p) Surface-tension-driven flow is not observed in water layers sub- 
jected to a propane flame exposure. Even though one would expect these 
flows theoretically, the inevitable contaminants of water precludes them 
in practice. 



VI. Conceptual Design of the Spacelab Experiment 


Description of the Experiment ; The objectives of pool burning experiments 
in the Spacelab are to make: (i) qualitative observations; (i) quan- 

titative testing of various gravity-related hypotheses; and (iii) quan- 
titative examination of the various assumptions made in theoretical 
studies. The experiments have to be such as to be uniquely possible only 
in the essential environment of the Space. Both the qualitative and 
quantitative studies will relate to flame spread and ignition aspects of 
pools of fuels at temperature lower than their flash point. These studies 
must lead to a distinct advance in the state-of-the-art and science of 
combustion of liquid fuels. 

Whereas the quantitative data are expected to enable us to verify the 
expected or speculated phenomena, the qualitative observations will serve 
two specific purposes. First, they aid to consolidate the quantitative 
measurements in some detail, say via photographic visualization. Second, 
they will reveal some features, phenomena, and processes which have not 
been identified before in our project research. 

The experiments would also establish the thresholds of validity of 
various concepts and notions, and existence of different limiting modes 
of flame spread and ignition. This would in effect lead to identification 
of limiting combinations of conditions which would permit a particular 
phenomenon to occur. Not only new useful information but also novel 
unfamiliar notions are expected to evolve out of these threshold 
delineations. 

The experiment to be carried out in the Spacelab must fulfill certain 
inherent constraints. The experimental set-up has to be large enough to 
permit the planned observations and measurements, and yet small enough to 
(a) permit visualization through limited -view-angle photographic means 
and (b) assure safety of the Spacelab and its crew in the event of a 
spill. The set-up has to be relatively insensitive to extraneous distur- 
bances (such as g-jitter). It should require little or no calibrational 
effort in flight; and must be as much automated as possible to avoid 
potential human error. 

Figure 10 shows a schematic of the conceptual Spacelab experiment. 

The heart of this experiment is the pan which holds the fuel. This tray 
is expected to be made of Pyrex. The dimensions of this pan are 2.54 cm 
wide, 2.54 cm deep and 35 cm long. The false bottom, whose height 
location can be adjusted by two thumb screw posts, is provided so that 
the test depth of the fuel layer can be adjusted to any value between 
about a few millimeters to about 2.54 cm. The idea of an adjustable 
false bottom is conceptually quite attractive; but in practice it may be 
much easier to replace pans of different depths. 

In the back of the tray, a ruler scale marked in mm divisions is 
attached. Sliding horizontally under the ruler is a shutter, made of a 
fireproof material such as asbestos, which when in position will snug gly 
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cover the 2.54 cm wide, 35 cm long tray. Upon completion of an experiment, 
this shutter is mechanically slid forward to extinguish the flame. 

Supply fuel is introduced into the tray through a port at the 
bottom. The false bottom will be so dimensioned as to loosely fit in the 
2.54 x3.54 cm tray such that the fuel introduced under it will slowly fill 
the entire tray to the desired free-board (of zero or over 6 mm) height. 

This scheme is expected to produce the least disturbances in the pool 
during the fuel-loading process'. The fuel itself will be carried aboard 
the Lab in containers (perhaps spherical) of inside diameter = 8 cm. This 
would provide enough volume (» 2/0 ml) of fuel to fill the tray 100% full. 

The containers have to be provided with an internal elastic membrane - 
introduction of compressed air to one side of which would expel the fuel 
on its other side into the fuel tray. Upon completion of the test, the 
used fuel can be drawn back into the container by permitting leakage of 
the air. Perhaps a well-designed cylindrical syringe is an alternative 
to this spherical shell concept. On one end of the tray, a pipette will 
be arranged in such a way as to gently deliver a few drops of an 
accelerant (such as methanol) whose flash point is low enough to be 
ignited by an electrical spark. The accelerant has to possess, addition- 
ally, a high dynamic viscosity to preclude rapid spreading over the test * 
fuel surface. It also should be immiscible in the test fuel. 

About midway along the length of the pool, three fine wire thermo- 
couples will be placed in the liquid, one at the surface, one at mid- 
depth and the third at the false bottom. A single thermocouple will also 
be used to continuously monitor the test chamber gas phase temperature 
before, during, and after the flame spread test. 

A flow tracer, perhaps finely ground particles of magnesia or alumina 
will be predispensed into the fuel supply container which can be shaken 
thoroughly before expelling the fuel into the tray. Besides the thermo- 
couple outputs, the majority of data of flow patterns in the liquid and 
gas phases and the flame spread itself are expected to be obtained on 
’film' or 'tape' via high speed photography. 

As far as study of preignition motion is concerned, an electrically 
heated rod of a known diameter is expected to be used as heat source. 

Its longitudinal location will be at the midpoint of the tray. The 
separati.on distance between the fuel surface and the ignitor is adjust- 
able. Motion picture records will yield both the time to ignition and 
the flow patterns preceding this event. 

3 

Procedure : The following major steps constitute the procedure of the 

flame spread experiment. 

1. Place the tray in the chamber. 

2. Adjust the false bottom to the desired depth. 

3. Install the thermocouple rake in the tray, and the gas phase 
thermocouple. 

4. Connect the fuel container to the tray. 
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5. Connect air supply to the container. 

6. Actuate the shutter to check its operation. Deactuate. 

7. Position the ignition rod and accelerant pipette in place. 
Actuate spark to check. 

8. Agitate the fuel container to mix the flow tracers. 

9. Admit air to fuel container to slowly introduce the fuel into 
tray to the rim. 

10. Turn on the thermocouple data acquisition system. 

11. Note the initial temperature of fuel, and chamber pressure 
and temperature. 

12. Turn the camera on. 

13. Note the settling of the fluid disturbances. 

14. At time-2, introduce accelerant. 

15. At time 0, strike the ignition spark. 

16. Observe and record. 

17. As the flame spreads all along the tray, activate the 
shutter extinguisher. 

18. Turn off the data acquisition systems and cameras. 

19. Release air from the fuel container to retrieve the spent fuel. 

20. Disconnect the container from the tray and the air supply 
from the container. 

21. Mark and stow away the spent fuel container. 

22. Vent the chamber. 

23. Go to step 2 above. 

'The preignition study experiment involves similar steps as follows: ‘ 


1-6. Same as above. 

7. Position the ignition source at the midlength in the middle 
of the tray rim, 

8-13. Same as above. 

14. Initiate energy flow to the ignitor. 

15. Set time = 0. 

16. Observe and record. 

17. When ignition occurs activate the shutter extinguisher. 
18-23, Same as above. 
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APPENDIX A 


Simulation Techniques 


There are five bench-testing methods useful to simulate gravitational 
effects. These are to: (i) use bodies of progressively smaller physical 

dimension [14-18]; (ii) use other extraneously applied and controllable 
body forces; (iii) use a centrifuge upon the arm of which the experiment 
is allowed, or made, to slide; (iv) use variable pressure in the fluid 
to alter gas-phase density [19-22]; and (v) use dyes in a water tank, 
the dyes being chosen variously buoyant relative to the working fluid, 
namely the water [23-25], 

These bench methods of simulation share some common concepts. Because 
gravity-dependent forces usually appear in our considerations in the form 
pgd J of the density or density differentials, gravitational constant and 
the linear dimension cubed, it is conceivable to control in an experiment 
the relative intensity of these forces by adjusting one or more of these 
variables. It is also possible to adjust gravity forces relative to other 
forces by suitably manipulating the magnitude of these 'other forces'. 

For example, the gravitational force may be reduced relative to the 
capillary force by reducing the Bond number Bo = pgd;2/0. This may be done 
by reducing p, gravity g or dimension d or by increasing the surface- 
tension a [15]. Porous wicks and sponges to hold liquids in any desired 
geometrical form illustrate the role of small pore size d which is in 
essence same as that of a reduced gravity constant g. Similarly, gravita- 
tional forces may be reduced relative to inertial force by increasing the 
Froude number Fr = pv2/pgd which may be done not only by reducing g and/or 
d but also by increasing the velocity V. A reduction in the Grashof 
number Gr = gd^p^AT/p^, as an additional example, reduces the gravity 
force relative to inertial and viscous forces. This reduction again may 
be accomplished equally well by reducing gravity g, or by reducing the 
density p [19-21] and the linear dimension d [15]. 

There are, however, practical constraints in employing indirect 
methods of simulating low-g conditions. Among these constraints in em- 
ploying indirect methods of simulating low-g conditions are the following: 

(1) Probing Difficulties : It is obvious that increasing difficulty 
is encountered in placing various measurement devices in an experiment 
as the scale of the experiment gets progressively smaller. The distur- 
bances caused by the presence of the probes and their supports also become 
progressively and unacceptably severe. 

(2) Physical and Chemical Limitations : As the density is progress- 

ively reduced to bring about the same effects as brought by the reduced 
gravity, certain physical and chemical limitations are approached. 

Examples of the physical limitations may be seen by considering the role 


played by pressure as a basic thermodynamic property of the working fluid 
and as an essential factor in the definition of a continuum fluid. 

Examples of chemical limitations may be identified by considering the 
role played by pressure in determining the reaction kinetics and mechanisms 
no matter whether they are simple or complex. 

Heat transfer mechanisms, which are relevant not only to maintain 
combustion but also to suppress combustion, also suffer with the intrusion 
of slip-flow effects (Knudsen free molecule flow effects) as the gases be- 
come rarified and non-continuum. If reduced gravity conditions were 
sought to be obtained by reducing the density differentials, a whole set 
of novel heat transfer problems are encountered in the limit of vanishing 
driving potentials. These novel problems pertain to systems only slightly 
perturbed from their equilibrium states* 

(3) Short Life-Times ; An examination of the time scales of interest 
in thermal response of the condensed -phase, vaporization, diffusion and 
mixing in gas-phase and combustion of fuels indicates that small-scale 
experiments last far too briefly to permit the required measurements. The 
availability of suitable measuring tools with acceptable precision and 
accuracy as well as fast response is a matter of separate but related 
concern. 


Consider the condensed-phase hydrostatics and dynamics, for an 
example. Suppose one wishes to employ Bond -number modeling. Then to 
preserve the Bond number pd^g/c? invariant with a particular liquid, d^g 
is to be kept invariant i.e., d«g“3-/2. if the Bond number considered is 
much greater than unity, that is in gravity-dominant regime, with Froude 
number « 1 the characteristic time t^Cd/g)!/^ which with the Bond number 
model becomes t a d3/2. Similarly, if Bo is much smaller than unity and 
Reynolds number is of the order of unity, that is in capillary-dominant 
domain, t a (pd3/a)l-/2 which for a particular liquid reduces to t<*d3/2 
again. Clearly this 3/2-power relation is valid for all Bond numbers. 
Then, a reduction of length scale by a factor of (1/10) reflects in. a 
characteristic time reduction by a factor of (1/10)3*2 or about (1/30). 
This reduction can pose formidable measurement, accuracy and response 
problems . 

Similarly, where the gas-phase convective, processes are attempted 
to be modeled by preserving Grashof number gd3 gATp2/p2 > gad-3 and the 
response time is of the order t^p/ (pgdgAT) so that t^l/gd^d^. A l/10th 
reduction in the scale of the experiment then cuts the time scale by a 
factor of (l/lO)^, that is by l/100th. By considering modeling of other 
aspects, in a similar manner, the headaches associated with reduction in 
the linear scale of the experiment to simulate reduced gravity may be 
clearly seen through the more drastically reduced characteristic time. 

(4) Visual Observations : Visual and photographic observations, 

likewise, also become immensely difficult or impossible due to the 
necessity of minute ^dimensions to be viewed in brief time frames. Speed 
and resolution, two often conflicting parameters of photography, are not 
always available to record an experiment. 



Other simulation techniques employed include (a) conducting the 
experiment in a freely falling chamber with various counterweights to 
achieve different g- levels; and (b) conducting the experiment in an 
aircraft whose flight trajectories are so prescribed as to realize the 
desired g-variation. 

Free-fall simulation of low-g has been reported in the literature 
quite extensively. The two most .noteworthy groups working in this area 
are from Japan [26-28] and IJS NASA-Lewis Research Center [29-34]. The 
time available in a free-fall facility is of the order t“hQ*5/2.2 sec, h 
being the height of fall in meters. A height of 5 meters thus gives 
about one second of fall time. In order to double this time, the height 
has 1 to be increased to about 20 meters; and to triple it, nearly 45 meters. 
The air resistance to the falling capsule becomes a severe problem with 
large fall heights. A double-capsule system, in which the outer capsule 
acts as a drag shield, is employed to overcome this problem. Another 
technique used is to drop the capsule in a near vacuum as is done in the 
Lewis' 500 foot facility. 

Whereas a free-fall gives an essentially zero gravity (i.e., 
gravity of the order 10”? g), counterweights may be employed to attain 
various fractional gravities (i.e., g lying between XO^ - ioO g Q ) u One 
of the most important things to remember in free-fall experiments is that 
the characteristic time of the phenomenon under study has to be less than 
the fall time. This requirement is not always met, especially in experi- 
ments dealing with liquid fuel evaporation, vapor/air mixing and combustion. 
Furthermore, transiting the experimental conditions from 1-g to 0-g> the 
disturbances are not quickly dissipated. As a result, the true initial 
conditions of the zero-g experiment are not well defined. Thus the finite 
and short duration of the low gravity ends up to be one of the mortal 
shortcomings of the free-fall technique. 

Aircraft executing Keplerian flight trajectories are used to 
obtain longer durations of low and programmed gravity conditions. Kimzey 
[35-39] used this technique to study combustion of various polymers. 
Free-floating capsules are normally used to isolate the experiment from 
the air frame. The practical aspects of flying the airplane pose time 
limitations which permit only about 10 seconds of true zero-g for experi- 
mentation. Maneuvering the aircraft in such a manner as to avoid 
collisions between the experiment capsule and the airplane walls is not 
the easiest of the tasks for the pilot even if one is satisfied with this 
short true zero-g duration. 
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APPENDIX B 


SOME RELEVANT PROPERTIES OF SOME LIQUID FUELS 
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Formula 

Mol. wt 

P 

g/cm 3 

Ah c 

kcal/g 

Ah v 

cal/g 

Si 

cal/g°C 

C Pv 

lean 
limit 
X vol. 

BP @latm 

°C 

FP 

1GN 
SponL . 
°C 

Viscosity 
Cp«3 0CJ 

Sur f . Tension 
dyne/ cm 
(<§ °C> 

1 . 

Methanol 

ch 3 oh 

32 

0.800 

4,74 

263 

0.57 

0.40 

5.5 

64.5 

15.5 

470 

0.82(0) 

0.60(20) 

0.46(40) 

0.40(50) 

24.5(0) 

22 . 6 ( 20 ) 

2 . 

Ethanol 

c 2 h 5 oh 

46 

0.790 

6.40 

200 

0.56 

0.46 

3.5 

78,5 

21.7 

392 

1.77(0) 

1 . 20 ( 20 ) 

0.83(40) 

0.59(60) 

0.50(70) 

24.1(0) 

23.0(10)* 

22.7,5(20)* 

21.90(30)* 

3. 

Acetone 

(CH 3 ) 2 CO 

58 

0.791 

7.36 

125 

0.51 

0.34 

2.1 

56.2 

- 9.4 

538 

0.40(0) 

0.32(25) 

0.28(41) 

26.21(0) 

23.70(20) 

21.16(40) 

4. 

n-Hexanol 

(ch 3 )(ch 2 ) 5 oh 

102 

0.814 

8.70* 

- 

- 

- 

- 

158 

74 

- 

- 

- 

5. 

n-Butanol 

<ch 3 )(ch 2 ) 3 oh 

74 

0.810 

8.60* 




1.45 

117.5 

43.3 

343 

5.2(0) 

3.4(15) 

2.3(30) 

1.4(50) 

0.93(70) 

26.2(G) 

24.6(20) 

22.1(50) 

6 . 

Pentanol 

(ch 3 hch 2 ) 4 ch 

88 

0.811 

- 

- 

- 

- 

- 

137.3 

40.5 

343 

- 

- 

7. 

n-Decane 

(cr 3 )^ch 2 ) 8 

142 

0.730 

10.56 

36 

0.52 

0.40* 

0.70 

174.1 

43.9^ 

208 

0,92(20) 

- 

8 . 

n- Nonane 

(ch 3 ) 2 (ch 2 ) 7 

128 

0.718 

- 

- 

. - 

- 

0.74 

151 

31.1 s 

206 

0.71(20) 

- 

9. 

n-Dodecane 

(CH 3 ) 2 (CH 2 ) 10 

170 

0.749 

- 

- 

- 

- 

0.60 

216 

74 s 

204 

1.35(25) 

- 

10 . 

Dipentene 

C 10 H 16 

136 

0,840 

- 

- 

- 

- 

- 

178 

68.3 

- 

~ 

- 

11 . 

o-xylene 

C 8 H 10 

106 

0.897 

10.30 

80 

0.41 

0.40* 

1.0 

144 

23.9 

464 

1 . 1 ( 0 ) 

0.88(16) 

0.81(20) 

0.63(40) 

30.1(0) 

12 . 

Kerosene 

-- 

— 

0.81* 

10.30* 

70 

0.46* 

0.40* 

0.6 

250* 

49 s 

260* 

- 

- 

13. 

2 -n- 0 ctanol 

(CH 3 )2(CH 2 )5 

CHOH 

130 

0.822 

9.70 

- 

> 

~ 

- 

86 

82.2 

- 

10.6(15) 

27.53(20) 

14. 

Benzaldehyde 

c 6 h 5 cho 

106 

1.042 

7.92 

- 

- 

- 

1.4 

178 

73.9 

192 

1.39(25) 

40.04(20) 

15. Benzyl alcohol C 6 H 5 CH 3 O 

Source: various handbooks 

* our estimate 
§ closed cup 

108 

1.042 

10.16* 

84 

0.38 

0.40* 

1.3* 

205 

104.4 

436 

5.8(20) 

27.7(10)* 

23.5(20)* 

27.4(30)* 



APPENDIX C 


Flame Spread Theory 


Consider a fuel layer of undisturbed depth h TO at an initial temperature 
situated in a chamber of pressure P a . Let this fuel's flash point be 
T 0 , density p, dynamic viscosity y, surface tension (with ambient air) a, 
thermal diffusivity a, and Prandtl number Pr (Figure 11). 

Suppose a flame is initiated over this liquid surface at x = - 00 and 
allowed to steadily propagate at a velocity U in the positive x-direction. 

We will measure y from the fixed bottom surface of the pool y=0, to the 
free surface y=h. x=0 is fixed on the steadily propagating flame. The 
fuel layer is of infinite width in z-direction. The x- and y-directional 
velocity components are u and v respectively, as shown in Fig. 11. In the 
chosen coordinate system, the flame front remains stationary while the 
fuel layer is fed backwards into the flame at a steady velocity. As a 
result of nonuniform surface temperature distribution near x=0, surface- 
tension gradients are invoked to cause perturbation of the flow from the 
uniform -U flow. 

The objective of our study is to obtain the flow and thermal fields in 
the fuel layer as dependent upon the various properties mentioned above, 
especially the flame spread velocity, surface-tension gradient, gravitational 
constant and flash point. 

In the following, we employ Boussineq approximation, i.e., density of 
the liquid is constant in all respects except in producing buoyancy. The 
steady state continuity equation is written as 


3x 3y 

Integrating from y = 0 to h, where h is the fuel depth, 

rh g 

1 dy + v (x) - v =0 
I 8x 3 s 7 o 


Since the base is impervious v = 0 so that 


3u . . ' ’ _ 

dy + v g (x) = 0 
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FIGURE 11: FLAME SPREAD SCHEMATIC 
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Tlie Leibnitz rule says 


■b(3) 

a(3) 


af 

83 


dy = 


33 


f (y.i 3)dy - f (b, 3) 


dt +£(a ’ B> ft 


so that Eq. (1) becomes 



udy - u 


dh 

dx 


u 


dh 
c 

o dx 


+ v = 0 
s 


( 2 ) 


The third term is identically zero since the base is rigid. We denote 
the x-component velocity at the free surface u(h) by u g . Now the velo- 
city in the y-direction at the free surface is composed of two parts: 
first, due to displacements in the free surface which under steady state 
becomes u s dh/dx by geometrical considerations and second, due to surface 
evaporation v ge . Neglecting evaporation, Eq. (2) becomes 

d 

— udy = 0 (3) 

dx 

o 

Equation (3) is valid irrespective of the variance (or constancy) of h 
with respect to x. Integrating Eq. (3), and applying the boundary con- 
dition that far ahead of the flame the velocity is uniform and equal to 

-u. 


■h 

udy = - Uh^ (4a) 


That is, the velocity distribution u(y) at any x will be such that the 
mass flow rate is independent of x. The momentum and energy equations are: 


u 



v 


3u 1 3P 

3y p 3x 




(5a) 



3v , 1 3P 
3y p 3y 


= 7 V 2 y - g + g3 (T-T M ) 

r 


(6a) 


u 


3T , 3T 
— + v — 
3x 3y 


aV 2 T 


(7a) 


All properties are considered invariant. Viscous dissipation of energy 
is ignored. is the Laplacian operator. P is pressure, y is dynamic 
viscosity of the liquid and a is its thermal diffusivity. In the 
following paragraphs arguments are presented to further simplify these 



conservation equations 


It will be assumed that the axial diffusion of vorticity and heat 
are both negligible. This requires that terms of the form 32 /g x 2 
neglected. In addition inertia terms in the momentum equations and 
convection of heat normal to the surface in the energy equation will both 
be neglected. Equation (5) thus becomes 




As the y-directional velocity component v and its derivatives are 
negligibly small, Eq. (6a) reduces to 


(5b) 


f- - ~ pg + Pge(T-Tj 


and Eq. (7a) to 
3T 


u 


3x 


= a 


3 2 T 

9y 2 


(6b) 


(7b) 


Equations (5b) - (7b) have to be solved along with Eq. (4a) in the 
domain 0£y<^h, 0<^x <°° under the following boundary conditions. 

At the base, y=0, since no slip is permitted, u = -U. We will assume 
that the base is a perfectly conducting boundary kept at a temperature 
(which is also the temperature of both the supply fuel and the 
atmospheric air). 


At the free surface y = h the pressure is assumed to be uniform and 
equal to the atmospheric pressure P a . The shear in the x-direction at 
the free surface has to be balanced by the surface-tensional stress, 
y3u/3y = do/dx, a being the surface tension of the liquid-atmospheric air 
interface. Furthermore, at y = h, we assume that the temperature gradient 
3T/3y is proportional to the difference between the local temperature T 
and the ambient temperature 

Also at the free surface but at x=0, the temperature of the fuel is 
equal to the flash point T c . Written mathematically these boundary 
conditions are as follows : 


y - o 

u = - U 

all x 

(8 a) 

y - 0 

T = Tc 

all x | 

(9a) 

y = h 

y3u/3y = da/dx 

all x 

(10a) 

y ~ h 

3T/3y = k(T-T 00 ) 

all x 

(11a) 

y = h 

►3 

II 

H 

O 

x - 0 

(12a) 

y - h 

P - Pa 

all x 

(13a) 
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All y 
Al l y 


u = - U 


x -> 00 


(14a) 


T - T m x-+“ 

The velocity u is nondimensionalized as u = ph u/y. The flame spread 
velocity U is specially denoted as Re = ph^U/y. The temperature is 
reckoned relative to the fuel’s initial temperature To, and normalized by 
the diff erence between the flash point T 0 and initial temperature. Too so that 
T = (T-Too)/ (Tq-Too) . x and y are measured in units of the pool depth so that 
x = x/hoo and y = y/h„,. The pressure is rendered nondimensional accord- 
ing to P = ph£ P/y. When one substitutes these nondimensional variables 
into the governing equations, the problem at hand reduces to the following: 


ap a 2 u 

85 3y 2 



(5) 

& = Gr(- 
3y 

1 __ 

- -i- + x) 
gAT J 


(6) 

3 2 T 

— 2 - Pr 
3y 

- 3T 
u — 

3x 


(7) 

•h 




udy = 

• 

O 

- Re 


(4) 

y - 0 

u = - Re 

all x 

(8) 

y - 0 

T = 0 

all x 

(<0 

y = 1 

3u/3y = ^3T/3x 

all x 

(10) 

y = 1 

3T/3y = kT 

all x 

(11) 

y ~ 1 

T = 1 

x - 0 

(12) 

y - 1 

P = P 

a 

all x 

(13) 

all y 

u = - Re 

X -V 00 

(14) 

all y 

T = 0 

x -> °o 
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In the above, several dimensionless parameters have evolved. These are 
defined in Table II. Specifically, Gr is Grashof number, the parameter 
which deals with the gravitational acceleration. 6 is the liquid's 
C specific volumetric expansion coefficient (1/°K) which when multiplied 
by AT = (Tq-T^) °K becomes gAT, nondimensional, a parameter specifying 
the expansion characteristics of the liquid* Pr is the liquid's Prandtl 
number, ratio of kinematic viscosity (y/p) to thermal diffusivity 
a(EK/pCp), Re is the flame spread Reynolds number. <p is the non- 
dimensional sensitivity of the surface-tension a to temperature T;_ 
o=da/dt, is assumed to be a given constant negative number; $sda/dT, iff 
being nondimensional surface-tension, k is simply (khco) where k is the 
proportionality constant in Eq. (11a). Finally P a is the nondimensional 
{ pressure of the atmosphere. 

It Is clear from the nondimensional_version of the problem statement 
that_the dependent variables u» T, and P are functions of x and y with 
Gr, <£, Re, Pr, gAT, k and P a as parameters, h, we stress again, is a 
function of x which we will here assume to be a constant. This is a 
highly questionable assumption, the seriousness of which can_only be re- 
vealed by experimental data from the Spacelab. In Table II h*l, i.e., 
h = h^ according to this assumption. 

Differentiating Eq. (6) with respect to x and Eq. (5) with respect 
to y, pressure can be eliminated; 

84 - Gr 4 (15) 

W 8x 

Equations (7) and (15) thus are the two equations for the two unknowns 
u and T. Suppose now that T is a linear function of x, an approximation 
which would imply that u is predominantly a function of y with its 
dependence on x being exhibited only through the constant 3T/9x at the 
free surface. If 

||=C(y), (16) 

9x 0 

substitution in Eq. (15) gives 

i 3_ 

^=Gr C 0 (y). (17) 

3y 

This means that the left hand side is a function only of y. The lower 
derivatives of u(y) as well as u itself are also assumed to depend on 
x similarly. Differentiating Eq. (7) with respect to x, 

3 3 2 T „ /- 3 2 T 9T 9u\ 
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TABLE II 


DEFINITIONS OF NONDIMENS IONAL QUANTITIES 


u = ph u/y 
00 

Re 5 ph u/y 
00 

a 5 ph ro cr/y 2 
= ph‘p // 

5 (T-TJ/CVV 

5 y/K 

s 4(T -1 ) ph / v 2 5 

' Q oo' ^ 00* n 

= kh 

00 

5 d ' /K 
00 

5 4'£4/pC p (T o -T m ) 

2 

t = yt/ph^ 


v = phv/ V 

I 5 shfp 2 /p 2 

op 

Fr = y/pa 

P 2 ph 2 P /y 2 

a ^ oo a H 

h - h/h 

00 

X 5 x/h„ 

♦ (V T - >h -/i* a Pr 

Gr = gh 2 p 2 6(T o -T„)/p 
A = A/h. 



Due to the preceding assumptions related to the linearity of T(x)» all 
terms in this equation except the first term, are zero. Hence 


2 — 
3 3T 


_2 

3y 3x 


= 0 


Integration with respect to y twice thus gives C (y) . 


C Q (y) = C ± y + C 2 


(18) 


where C-, and C» are integration constants. At the bottom of the pool, 
T is unxformly equal to zero for all x. Thus Eq. (9) requires that 
3T/3x be zero at y = 0, so that C 2 = 0. 


C Q (y)- C x y 


(19) 


Substituting Eq. (19) in Eq. (17) and integrating three times, we obtain 
the velocity. 


-4 -2 

" = Gr C 1 47 + C 3 2T + C 4 y + C 5 


( 20 ) 


C^, and C,_ are integration constants, constants as is a constant 
to the extent 3T/3x is independent of x. 

Integrating Eq. (19) with respect to x. 


T = C y x + f (y) 


( 21 ) 


where f (y) is a function to be determined from Eq. (4). Substituting 
Eqs. (19) and (21) in Eq. (7), 

d2f ~ ' 

72 = Pr c ! uy 

dy 

which can be integrated wth u(y) known from Eq. (20). The result with 
Eq. (21) gives T. 


T=Pr C 


i[ 5GrC i fr + 3 sIt + 2C 4 It + S fr] + < c 6 + c 1 ; )5 :+0 7 (22) 


The six constants C 1 ,C 3 ,...C 7 are now to be obtained from Eqs. (4), 
(8)-(12). Of special interest is that if Eq. (11) has to hold for all 
x, k is required to be unity. (This would correspond to a unit Biot 
number situation in thermal conduction.) 



40 Re 


1 

C 3 “ I (*'^ Gr ) C 1 

C 4 


= - Re 


Pr Re C, Pr C 


C, » 1 + 


C = 0 


' 1 / 23 0 -\ 

mo- (7Fo 0r * ♦) 


Applying these constants to Eqs. (20) and (22) and sorting out the terms 
containing the Grashof number and <j>. 


u = - Re + C 

T = (1 + Ci) 

Pr C 2 
+ ”240 


(£ - ± y 2 +j- V-i/2. y 2 -yVl 
8 \ 3 10 y 15 ' / 2\2 y Vj 


(23) 


[ 


— (100y 7 - 567y 5 + 490y 4 + 99y 3 - 122y) 


+ $ (9y~* - 10y 4 - 3y 3 + 4y) 


■ 

> 


(24) 


The pressure field can be found from an integration of either Eq. (5) 
or (6) and boundary condition of Eq. (13). We will decline from_this 
derivation and focus on the implications of the predicted u and T fields. 


Figure 3 shows the velocity profile given by Eq. (23) for three 
different values of the ratio Gr/<j>. This ratio is an indication of 
whether the gravity effects or the surface-tension effects dominate in 
the convection process. A pool of n-butanol of 0.02 m depth, supply 
temperature T^, = 299.7 K, will have a Gr of 6000. The surface-tension 
parameter $*-100. Thus the Gr/$ =-60 curve indicates the earth-based 
experiment while Gr/$ = 0 is the zero gravity experiment. Also shown in 
the figure is the profile when Gr/$ •> -®° . Several interesting points are 
evident from Eq. (23) and Fig. 3. 
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1. Far ahead of the flame where no surface-tension and temperature 
gradients exist the fluid velocity will be unform u = -Re. At smaller 
values of x, the velocity profile is perturbed by the cooperative action 
of buoyancy and surface-tension effects in such a way that at about y-1/3 
a velocity maximum in the negative x-direction would oCcur. With zero- 
gravity (i.e., with only surface-tension effect) this maximum occurs 
exactly at y = l/3 with u = - Re + C], $ /12. As the gravity effects be- 
come relatively larger, tfie depth at which this maximum occurs will slightly 
decrease towards y« 0.274 when Gr/4> -*■-<» where u~ -2,319 Re. 

At about y = 2/3 (as at y = 0), u = -Re. With increasing Grashof 
number, this value of y where u * -Re will decrease towards about 0.59 as 
Gr->°°. At the surface y = 1, u = u g = 2. 33 Re when Gr/<j> ->■ 0 and u^Ug* 1.121 
Re when Gr/<j> In general, 


°s°¥( 7i -ii Gr ) 

Another important feature of Fig. 3 is that As Gr/$ gradually 
increases towards the slope 9u/9y at the surface gradually decreases 
towards zero. This fact is quite evident from the boundary condition 
Equation (10) . 

2. In the present theory_, the coefficient is independent of x as 
a result of the assumption 3T/3x is a constant with respect to x. In an 
exact solution these quantities will be functions of x. Then there exists 
a certain x at which 



and where the local surface velocity u g is zero. If x is larger than this 
certain stagnation point, fluid elements at the surface move towards the 
flame (i.e., in the negative x-direction) whereas behind the stagnation 
point, the surface velocity will be positive. It is this forward movement 
of fluid induced jointly by the surface-tension and gravity effects which 
would convect heat forward to enable the flame to spread, 

3. The linearity of temperature with respect to x is not merely an 
ad hoc assumption. For unexplained physical reasons, this simple behavior 
has been clearly demonstrated by the detailed numerical solutions of 
Torrance [8,9]. This linearity makes the x-dependency of the velocity 
profile rather weak. As long as our region of interest lies within 
0£ x _< - 1/C-p the velocity profile retains a fixed shape; to the right 
of this domain, the velocity will he uniform and equal to -Re. The 
abruptness with which the linear temperature distribution would indicate 
the velocity field is a quirk of the present theory. This issue requires 



further study. 


4. Note from Eq. (23) that when the. gravitational effects are absent, 
the functional form of the velocity profile is precisely that which has 
been obtained before by such authors as Levich [40], Adler [42], Yih [41], 
and Sirignano and Classman [7]. L 

5. Figures 4 and 5 show the temperature profiles of Eq. (24) for two 
different flame propagation Reynolds numbers. Three different Grashof 
numbers are considered in these plots keeping the liquid Prandtl number ^ 
and surface-tension parameter constant. As stated earlier, a Gr of 6 x 10 
refers to earth-based experiment on a 0.02m deep pool of a fuel such as 
butanol at room temperature. The maximum temperature attainable in the 
liquid layer is T = 1 which occurs at the free surface y = 1 near the 
fjLame's leading edge x = 0. The lowest possible temperature is, of course, 

T = 0. 

Eventhough our theory considers the flow problem as fully developed 
the temperature profiles indicate a boundary layer type development. De- 
pending upon the conditions of Gr and Re, the thermal penetration occurs 
to a depth of^ 5 from the surface in such a way that for y 1-6 the 
temperature T is zero. The magnitude of 6 depends on the distance ahead 
of the flame front. For x ,x 0 (=- l/C^) the temperature profile is 
uniform and equal to T_= 0 (see Eq. ( 24 ) ). At x = x m ,6=0. As one moves 
towards x=0 from x = x^ gradually, 6 would gradually increase. 

An examination of Eq. (7) points out that if_u were a constant, 

(say u ), similarity solu tions must be possible for T(x,y) in such a way 
that S y/ /[(xoo-x)7Pr u s ] is the similarity parameter . In order to test 
this hypothesis, the calculations of Figs. 4 and 5 are examined to see 
how the depth of penetration would vary with x. The result is as follows: 



The constant is given in Table III for the different Gr and Re. 

Remarkable is that the constant is about 6 for the two larger Grashof 
numbers irrespective of the Reynolds number. It is about 2.5 for Gr=0. 
Eventhough one would wonder if the use of the surface velocity u g (which 
is a function of Gr, ^ and Re) in composing the similarity parameter is 
appropriate, the result is interesting. Further study of this aspect 
is required. Suffice it to note that under otherwise same conditions, the 
zero-gravity thermal field penetrates less deep into the fuel layer. 

6. At a low flame propagation speed, as shown in Figure 4, a Gr of 
1()5 results in a temperature field which is fully developed, i.e., the 
thermal penetration occurs all the way to the bottom of the pool. The 
profiles are nearly linear. Considerable heat losses may then occur to 
the base of the pool. At a medium Grashof number value (* 6xl0 J ) the 
penetration occurs to the base but only close to the flame front. However, 



TABLE III 


TEST OF BOUNDARY LAYER HYPOTHESIS 


Pr = 45 

<j> = -100 

< fsl "V 1/C 1> 

GR 

Re 


10 5 

2 

_ * 


10 

6.04 

6x10 3 

2 

5.94 


10 

6.05 

0 

2 

2.50 


10 

2.77 


* 

Fully developed thermal field. 

at zero gravity, the maximum depth to which the thermal perturbations 
penetrate is only about a quarter of the pool depth. 

7. A comparison of Figures 4 and 5 points out that while the general 
features of the thermal field remain changed, a faster propagating flame 
would result in a shallower penetration of the thermal disturbance. Even 
at the largest Gr considered, only the top half or so of the fuel layer 
will be influenced thermally. For lower Grashof numbers the temperature 
changes are limited to a thin layer of the fuel near the surface. 

8. The surface temperature variation evident in Figures 4 and 5 is 

specifically presented in Figure 6. The linearity arises from our basic 
assumption. From T s = l_at x = 0, the temperature decays to T g = 0 at a value 
of x = = - 1/Cy* - 3(<j>-llGr/140)/40 Re. Larger Reynolds numbers, smaller 

surface-tension parameters and smaller Grashof numbers result in a shorter 
region of thermal perturbance. Especially the gravity level decreases the 
extent of the x-region of activity by over an order of magnitude. This is 
an important result which would merit verification by experiments in the 
low-gravity. Figures 7 and 8 illustrate the variation of with Gr, Re 
and <j> . 

9. The extent to which this theory is compared with experiments 
involves a comparison of the velocity profile and measurements of flame 
spread rate as dependent on Gr and <j>. These are described in Appendix E. 

The measurements show that typically for alcohols, <j> ~ - 100 to 200, 



Pr s 45-100, Re »l-20, with Gr s 6000. The results discussed above 
lead us to conclude that a 0.02 m fuel layer will experience thermal 
perturbation to a depth of 0.005 m at zero-gravity. The length then 
is less than about 5 so that a physical length of 0.1-0.15 m is required 
to observe the region of thermal influence. Thus we can conclude that a 
tray of depth 0.02 m and length 0.4mwill suffice. As the experiments 
of Appendix E show, this length would permit exclusion of about 0.07 - 
0.15 m for the ignition transient decay and 0.05-0.10 m for the 
alleviation of the upstream boundary effects. 



APPENDIX D 


Ignition of Fuel Pools that are 
Initially Below Their Flash Point 


Consider a fuel pool whose initial temperature is below its flash 
point. Suppose an electrically heatbd wire, to be considered as a cylinder, 
is placed parallel to the pool surface at a distance D from the surface in 
the gas phase. Let the power (q* per unit length) be switched on at time 
zero. A transient free convection flow is thus initiated in the gas-phase 
surrounding the cylinder. A temperature field would develop to steady state 
in a short time. As the flow develops and thermal wave penetrates, a time 
will come beyond, which the liquid surface gets into thermal contact with 
the heat source. Suppose that the gas phase thermal field around the 
ignition source is cylindrically symmetrical (as expected in the absence of 
gravity). Suppose then that the isotherm of significant influence on the 
surface is of diameter d, this isotherm cutting the quiescent liquid surface 
resulting in a segment of length d T which is the width of heated liquid 
surface strip, d’ depends on D and d according to 




As a result of this heating, a surface temperature gradient is established 
as one marches along the liquid surface normal to the ignitor. The conse- 
quent surface-tension gradient would initiate a convection process transiently 
forming a pair of symmetrically placed cylindrical flow cells. If q’ were 
dispersed in the gas phase around the ignitor symmetrically, the fraction 
of q’ delivered to the liquid is given by 


f = 2 ~ sin' 1 (d’/d) 

Given this situation, to determine the flow patterns, the heat transfer 
distribution, the resultant production of fuel vapors, their dispersion in 
the gas-phase and ignition of the mixture so formed consitute an important 
problem both in science and practice. 

Our considerable search into the engineering, physics and chemistry 
literatures pointed out only a few references to this problem 12,4] all 
from Princeton. We therefore develop in the following a brief analysis 
to identify the process components of the problem. 

The quantity of energy f V delivered to the liquid is transported away 
by convection. The fluid velocities both in the x and y directions, u and 
v, are of the same order of magnitude in this ignition problem. The con- 
vective transfer of energy hence is described by 

fq’ - d’ pu C p (T-TJ (1) 

where d’ is the area of incident heat per unit width of the pool, pu ~ pv 


mix 





"wsyris 
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is the liquid mass flow rate, C is specific heat of the liquid of initial 

p 

temperature T^. T is the time- and space-dependent temperature of the 
liquid fuel. The velocity u is dependent on time and space as well. 
Assuming that the flow and thermal fields react instantaneously to the 
Influence of the ignitor-heating, we assume that u • u g , the 
x-directional velocity at the free surface of the fuel. From the flame 
spread theory (Appendix C, Eq. (23)), 


u 

s 


"fl /Gr 
4 y20 



( 2 ) 


Here $ is the sensitivity of surface tension to the fuel temperature 
do/dT. Implicit in Eq. (2) is that the surface temperature is representa- 
tive of the convective dispersion intensity of energy. Substituting Eq. 
(2) in Eq. (1), 


f<i’ =d'pC p u(T-T J * d’pC p u s (T s -T oo ) (3) 

which becomes 



with the nondimen sional quantities defined in Table II. Integration with 
T = 0 at x = A gives 



The x-directional extent of thermal penetration A is an indication of the 
cell dimension, ah yet undetermined function of time. Before further 
dealing with this quantity, the temperature averaged over x * 0 to A is 
obtained from Eq. (4) to be 



As the energy is thus convected, it is employed to heat two approxi- 
mately cylindrical masses of liquid of diameter A and of temperature 

?? x 


av 


fq' 


2 . 


ttA 

4 


dT 


PC. 


av 


P dt 


( 6 ) 


*1' 1 


54 



The factor 2 on the right hand side is to account for the two symmetrical 
cells. ttA 2 /4 is the volume of the cell per unit width of the pool. With 
the nondimensional time t defined in Table II, Eq. (6) becomes 

dT 

q « 2uA 2 — 321 (7) 

dt 

From Eq. (5) the temperature- time derivative is 


av 1 f 2 


1 pi_f /2 s - 1/2 M 
3 I A l— - *\l dt 


with which, Eq. (7) may be rewritten as 
r or nl/2 3/2 >7 


Mm 


Upon integration with the initial condition that there exists no cell 
(A * 0) at t = 0, 


= foT 

4tt 
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Equation (4) in conjunction with Eq;‘ (8) gives the surface temperature 
as a function of x. Equation (5) gives the average temperature of the 
cell to depend on time as 
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From Eqs. (2) and (4) the surface velocity is 
r r/Gr rVil/2 


2d (A-x) 


Equations (4), (8), (9) and (10) give a detailed structure of the 
transient flows developed in the ignition process. Ignition is expected 
to occur when the temperature a_t x=0 reaches the fuel's flash point 
temperature T 0 . Thus setting T=i and x * 0 in Eq. (4), with A given by 
Eq. (8), the ignition time t is obtained. 
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The extent of the cell at this time is 



( 11 ) 


( 12 ) 


Since for typical alcohols <J> lies in the range -100 to -200, q z 1.5-15 

(for q' * 10 _1 -10° kw/m, f = 0.5), d * 0.04 - 0.20, A* will be about 
0.33-15 cm, t* * 0.125 -10^ sec. at zero gravity and A* 2 10-15 cm and 
t* * 1-103 sec at normal gravity. 

In Figure 9, the predicted dependency of t* and A* on Grashof number 
and surface-tension parameters are shown especially attempting to point out 
the influence of reducing gravity to zero. Notable qualitative agreement 
is obtained between this theory and the experimentally observed flow 
patterns. Implicit in both the flame spread theory of Appendix C and 
Ignition theory of this appendix, we must stress, is the assumption of 
negligible hydrostatic depression of the liquid layer free surface. The 
impact of this assumption can be best assessed by comparison of normal 
gravity data and results obtained by experimentation in Spacelab. 



APPENDIX E 


Experimental Work 

The experimental part of our project was intended to check some of 
the theoretical predictions related to the feasibility and justification 
issues. These issues pertain to the dimensions of the pool, the choice 
of fuel, possibility of using inert liquids to examine the physics of flow, 
the choice of flow tracers, etc. 

1. Preliminary Experiment : Our preliminary experiments involved a 

simple cast iron tray of constant width and one meter length. Washed 
sand was placed in the tray to produce a flat inclined bottom such, that 
when the tray is filled with liquid fuel, the pool depth varied linearly 
along the length of the tray with zero depth on one end and 3.2 cm on the 
other. Some interesting observations were made with this apparatus. 

(i) When ignition was attempted by bringing a match flame to 
the vicinity of the liquid surface, ignition was found to be extremely 
difficult in the deep end. This is because of the excessive energy 
dissipation into the liquid phase as a result of strong circulation of 
the fluid produced by the surface-tension effects. Ignition in the 
shallow end, however, would occur readily, perhaps because of the 
suppressed liquid motion. 

(ii) Once ignition is accomplished in the deep end, even if with 
some difficulty, the flame would propagate towards the shallower end, 
first steadily and then in a somewhat faltering way. A flame initiated 
easily on the shallower end, however, would propagate but infinitesimally 
slowly; for all practical purposes, the flame burns there over a wick 
formed of the sand. 

(iii) If ignition is attempted somewhere in the middle of the 
pool length, the deeper the pool the more difficult the ignition and the 
easier the flame spread. 

These simple experiments pointed out the importance of surface- 
tension-driven liquid convection which plays opposing roles in ignition 
and flame spread. In ignition, it disperses energy away thus delaying 
the arrival of the instant when the liquid would produce enough vapors 
conducive to ignition. In spread, the convection offers the controlling 
forward heat transfer mechanism. A critical depth of about 2.5 cm was 
discovered in our preliminary experiments to permit formation of a 
convection cell which is unhampered by the lower solid boundary of the 
pool. This observation appeared to be independent of the particular fuel 
under test. 

2. Choice of Fuel : The tested fuels include kerosene, xylene, 

ethanol and butanol. Among these, kerosene is the only fuel which is a 
mixture of several basic hydrocarbon fuels. Even though the various other 



fuels have a stated purity (ex: 95% pure, tech, grade, or reagent grade, 
etc.) the impurity in most cases is caused by the presence of different 
isomers of the same molecular formula. Our constraints in the choice 
of the fuel included: easy availability, safety in view of toxicity, 

cleanliness in terms of sootless burning, lack of corrosivity, etc, 

Choice of a fuel which burns without soot is also desirable in view of 
reducing the effects of radiant heat transfer in flame spreading. And, 
of course, the fuel has to possess a relatively strong dependence of 
surface-tension on temperature. The flash point temperature has to be 
larger than the room or fuel supply temperature to accentuate the role 
played by surface-tension but not so large as to yield a pulsating pre~ 
cursor flame in the flame spread experiment . 

Our experiments rapidly pointed out that xylene (and other fuels 
with cyclic structure) are undesirably corrosive to the. pool pan and the 
associated sealants. Kerosene is too commercial a mixture to have well- 
documented properties. Ethanol and some other fuels are too close in 
flash point to the fuel supply temperature. Decane is too expensive to 
permit repetitive experimentation. With these observations, n-butanol has 
been chosen as our test fuel. 

3. Tray construction : Figure 12 shows the detailed construction of 

the fuel tray for the majority of experiments. Trays are fabricated of 
0.080 gage aluminum in such a way that the side walls are only l/2 M high. 
1/8" thick plate glass walls are then inserted in place with silicone 
sealant applied between the glass and aluminum to avoid fuel leaks. In 
some tests, the silicone sealant was replaced by a folded teflon tape 
and the glass walls are press-fitted. Xylene fuel attacks the silicone 
sealant seriously and the teflon tape less seriously. Six different 
widths ranging from 1/2" to 4" are thus accomplished for the finished 
pool trays. All the trays are 48" long except the smallest one which is 
24" long. All are 2" deep. The end walls of all trays are unlined 
aluminum. Each tray is provided with fitting false bottoms which are also 
fabricated out of aluminum; these are to provide a variable depth to the 
pools while at the same time resulting in a solid bottom . 

MacKinven et al [6] employed a substrate layer of water over which 
the fuel is situated in the desired thickness. The solid bottom arrange- 
ment chosen in our experiments is expected to eliminate difficulties of 
stipulating the velocity and thermal boundary conditions at the base in 
the theoretical analysis. 

A screw-hole was tapped in the bottom of each tray to accommodate a 
machine screw, removal of which after a test facilitates emptying the 
pan into a spent-fuel reservoir. A 15 cm ruler was situated in the back- 
ground for Obtaining the temporal location of the spreading flame front 
from video records. This entire tray set up was placed inside a 4* x4' 
by 8' long test chamber to avoid spurious drafts and to keep under control 
the combustion products. Various lighting distributions and backgrounds 
are employed to obtain the best records of fuel convection patterns 
photographically and via video taping. 




1. Aluminum tray frame 

2. Side glass 

3. End glass 

4. False bottom 

5. Silicone seal 

6. View line 


FIGURE 12: TRAY CONSTRUCTION DETAIL. 
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In a few experiments, a copper tube coil, carrying water at a chosen 
temperature, is placed in the space under the false bottom in order to 
achieve a prescribed initial temperature of the fuel. Due to difficulties 
associated with achieving a steady uniform temperature of the fuel by thus 
cooling or heating, the scheme is aborted in later experiments. 

4. Procedure . Typically, upon setting up the tray of certain width 
in place, the false bottom is installed to obtain the desired fuel depth. 

The fuel is introduced slowly until the free surface is either flush with 
the lip of the tray or below the lip to the desired free-board. Sufficient 
time is then allowed to permit the fuel layer to achieve quiescence during 
which time the tracer is prepared and lighting and background adjusted for 
photography. 

a. Flow Visualization : Several tracers are tested for flow 

visualization [13,44] : these include dyes (milk, food coloring, ink ? 

methylene blue, and others) and particles (dust, crushed aluminum shavings, 
talcum powder, crush'd polystyrene, fine powder of asbestos, and others). 

The asbestos powder when dispersed in the liquid yielded an unexpected 
effect. Fine bubbles of air desorbed by the asbestos rise in the liquid 

in the form of a string. This string of bubbles is deflected by any 
disturbing currents in the liquid, however mild. Whereas a string of 
bubbles is visible to the naked eye in the shape of the local instantaneous 
velocity profile, it is extremely difficult to capture this profile on 
film or video tape. 

Particles which float over the fluid surface are excellent to observe 
the surface flow patterns. But they may alter the surface-tension 
phenomena considerably. The most desirable tracer is one which is neut- 
rally buoyant in the liquid and yet is visible to the eye and camera. We 
have discovered that finely ground chalk or talcum of the color of one's 
choice (purple and blue appeared to us desirable) when soaked overnight 
in the fuel of concern (in such proportions as about 5% talcum to 95% 
fuel) would form a fine paste whose properties are expected to be similar 
to those of the fuel itself. This slurry could be applied to a quiescent 
pool of fuel with an eye-dropper specially equipped with a fine wire to 
convey the slurry to the liquid with least resultant inertial disturbance. 

This is the tracer technique we have chosen to employ in most of our 
experiments. 

b. Ignition Methods in the Flame Spread Experiment ; Three different 
methods of ignition are employed in the flame spread experiments. The 
simplest and the most extensively used one involves bringing a kitchen 

match flame close to the fuel surface (» 0.5 cm separation) and holding 
it there until ignition would occur. This method has the merit of pro- 
ducing the least amount of disturbance in the gas and liquid phases. But 
it also has the shortcoming of being ineffective when the flashpoint of 
the fuel is much higher than the supply temperature. To handle this 
situation, the second method is developed. In this, a small quantity of 
low flash point fuel such as methanol is gently relea.sed on to the test 
fuel surface at one end of the tray to form a thin floating layer which 
can be ignited easily with a match flame. If the quantity of accelerant 



released is large enough, the test fuel ignition occurs quite smoothly 
and a spreading flame evolves. In the third method, which is the least 
used in our experiments and which is meant for situations requiring con- 
siderable quantity of the accelerant because of large differences between 
the flash point and supply temperatures, a solid barrier flag fabricated 
out of masonite is employed to contain the relatively massive accelerant 
layer from spreading over the test fuel surface. Upon ignition of the 
accelerant, the barrier is gently removed. The liquid is expected to be 
disturbed to some extent during this removal, no matter how gently it is 
carried out. Whereas the mode of ignition did not cause any significant 
differences in the flame spread measurements in the 48" long trays, it is 
expected that the first method is the least disturbing and the third 
method is the most disturbing to the condensed phase quiescence. 
Correspondingly, the portion of the tray length influenced by the transience 
of ignition is expected to be largest in the third method. 

c * Heat Source in the Ignition Experiments : Two different heat 

sources are employed in the study of preignition flow patterns in the 
liquid. A heated nichrome wire of diameter 0.08 mm (B&S gage 40, 
resistance c 2.3 ohms/cm) is placed in the direction of width between 
ceramic terminals. The distance between the wire and liquid surface can 
be accurately adjusted. The advantage of this method is that the power 
input can be prescribed quite accurately. The disadvantage lies in that 
for any significant thermal communication between the relatively thin 
heated wire and the liquid surface, the wire has to be brought quite 
close to the surface, so close that even minor disturbances and vibrations 
in the building would cause the liquid fuel to touch the wire whereupon 
the surface-tension would keep the wire in continuous contact with the 
liquid . 

To overcome this difficulty, a propane gas diffusion flame fed 
by a feeble flow of the gas through a 1/16" dia. copper tube is devised. 

The resulting flame is blue; its size can be adjusted from about 5 mm 
to 10 mm, with its length about 5 times the diameter; it can be brought 
close to the fuel surface; there exists a thermal communication between 
the nearly spherical bottom-half of the flame and the fuel surface even 
at such large separation distances as about 10 mm. The power input rate 
and flame size are controllable by adjusting the fuel flow rate but some 
ambiguities arise due to the incomplete combustion, a severe characteristic 
of all diffusion flames. In most of our pre-ignition flow studies, the 
gas flame ignitor was employed. 

5. Results and Discussion ; In the following, the results of experi- 
ments are presented first on flame spread and then on ignition and their 
implications are examined in our concern of the feasibility issues. 

(a) General Observations : With all fuels the flame spread is 

smooth only in a certain range of supply temperature. For butanol this 
range is between 27.7 °C and 37.7 °C. If the supply temperature is lower 
than about 27 °C, the flame would propagate with pulsations, its leading 
edge propagating ahead to 2. or 3 centimeters and then withdrawing back. 
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To measure the flame speed consistently is nearly impossible when these 
pulsations occur. If the fuel temperature is greater than about 37 °C, 
the flame propagation is extremely fast (= 50 cm/s). Meaningful obser- 
vations of the liquid and gas phase processes are again extremely 
difficult. 

In the intermediate range of supply temperature (27.7 to 37.7 °C) 
the flame would propagate in an orderly and uniform fashion permitting 
careful observations. Normally the flame fronts in 4" wide pools are 
nearly flat with parabolic shape as the sidewalls are approached. The 
liquid phase convection is clearly visible both by bubble tracers and 
talcum paste tracer. Typical velocity profiles in the liquid phase ob- 
tained by this later technique are consistent with what is expected from 
the hypothesis of Appendix G related to the induction of a surface-tension- 
driven flow. 

In study of the motion of liquid in preignition heating, flow 
fields are obtained in a similar way. Because of the irregularities pro- 
duced by the surface-tensional attachment process between the wire and 
the liquid surface, we have aborted the idea of using heated wire as 
ignition source and resorted to employing a 5 mm dia. propane diffusion 
flame for the purpose. Depending upon the separation distance between 
the blue flame (bottom part) and the liquid surface, a transient convec- 
tion cell is developed. 

(b) Influence of the Tray Construction Details : MacKinven [6] 
pointed out that the glass liners along the sidewalls of the tray resulted 
in a consistently larger flame spread velocities compared to the velo- 
cities obtained in unlined aluminum trays. This is expected to be due to 
heat losses to the side walls from the propagating flame and convected 
fluid. We expect, as evident from the remarks below in relation to the 
pool width effects, that these heat loss effects will have to be taken 
into account in a comprehensive theoretical program to develop a correc- 
tion factor which would enable inference of flame spread in an infinitely 
wide tray from spread rates in trays of finite width. 

That the end walls are bare aluminum surfaces is deemed irrelevant 
so far as our ignition and flame spread studies are concerned. 

(c) Room Conditions : Most of the experiments reported here are 

conducted during spring and summer of 1978. The atmospheric pressure and 
humidity appear to be of no measurable consequence on the flame spread 
rate. The atmospheric temperature, however, appears to influence the 
spread rate via the fuel supply temperature. The reason for this is that 
the fuel is stored in a cabinet in which the temperature is essentially 
the same as that of the room. 

(d) Fuel Ty e : The nature of fuel effects the flame spread 
velocity strongly at a fixed supply temperature in a given pan. At an 
initial temperature of 20.5 °C, in the 5.08 cm wide 90 cm long pool with 



a fuel depth of 2.54 cm, the flame spread rates for xylene, kerosene and 
butanol are respectively 2.8, 0.6 and 1.2 em/s respectively. The flash 
point temperatures (open cup) of these fuels are respectively 23.8, 49, 
and 43 °C. Mixtures of kerosene and xylene are also tested in the same 
tray keeping the initial temperature 21 °C + 3%. The measurements are 
presented in Figure 13. The decrease in flame spread with increase in 
kerosene content of the mixture is not quite linear. 

Another set of data obtained on butanol has the objective of 
determining the influence of fuel purity on flame spread rate. Figure 14 
shows the spread rate as a function of the supply temperature for fresh 
butanol and for previously used butanol. There are obviously some 
significant differences, especially at lower temperatures . The spent 
fuel flame spread rate is, in general, lower than the fresh fuel flame 
spread rate. The video records lead us to believe that the surface- 
tension-driven currents in the liquid phase are not as intense in the 
spent fuel as in the fresh fuel. Ignition by match flame appeared to be 
a bit easier with spent fuel than with; fresh fuel. 

These observations lead us to believe that the spent fuel sur- 
face tension o and its sensitivity to temperature do/dT are lower than 
the same for fresh fuel. During combustion, the fuel surface is expected 
to encounter such products of combustion as soot and water vapor which 
are readily absorbed by the fuel. Knowing that even minor quantities of 
impurities in the fuel would enormously change the surface-tension 
characteristics, the observed lower flame spread rate and easier ignition 
of spent fuel are both logically expectable. 

Finally, it is interesting to note that xylene is found in our 
experiments to be rather undesirably corrosive, even to silicone sealant, 
to such an extent that an uncontrolled spill fire hazard had arisen. 

(e) Fuel Supply Temperature Effect : Figures 15,16 and 17 

present the influence of the fuel's initial temperature on the flame 
spread rate in trays of different widths and different fuel layer depths. 
Evident from these figures is that the flame spread rate increases from 
about 1 cm/sec at low temperatures to over 100 cm/sec at temperatures 
approaching the flash point (43.3 °C). The sensitivity of spread velo- 
city to variations in fuel's initial temperature is much stronger at 
higher temperatures than at lower temperatures. This relationship appears 
to be nearly exponential. For fuel temperatures above 38 °C or closer to 
the flash point, the flame propagation is so fast that even the motion 
picture records yield data with a great deal of scatter. 

(f) Free-Board Height Effect : The experiments reported here 

are all conducted with the fuel filling the tray nominally to 100% full; 
i.e., the fuel surface is nominally at the top edge of the tray walls. 

The nominal qualification is made to recognize that the fuel surface may 
be curved due to the surface-tension wetting of the walls. Nevertheless, 
the trays are so filled as to maximize the fuel loading without permitting 
spillage over the edges both before and during the flame spread. 
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The issue of how full is full is relevant in the viewpoint of 
not only the fire safety in the conduct of this experiment but also the 
energy feed forward mechanism itself. If the tray is only partially full, 
heating of the exposed part of the side walls, (the height of this exposed 
part being called 'free-board'), is expected to alter the mechanism of 
heating the fuel. Furthermore, the flame propagation in a partially 
filled tray would essentially resemble propagation in a narrow channel. 

The side walls then are expected to restrict the arrival of oxygen to 
the vapor and thus alter the flame spread characteristics. 

The net effect of the heat transfer and air flow patterns 
altered by the free-board walls according to MacKinven's experiments [6], 
is to reduce the flame spread velocity. As the free-board height is 
increased from 0 (i.e. , fuel surface flush with the top of the trays 
walls) to about 6 mm, the flame spread rate falls sharply. Further in- 
creases in free-board does not significantly influence the spread rate, 

(g) Effect of Tray Width : The data of Figures 15 to 17 also 

indicate the effect of tray width on flame spread rate. Figure 18 is a 
cross-plot of these data at one particular depth of the pool, viz: 2.54 
cm. The points in this figure are not necessarily actual data points 
but readings from Figures 15 to 17. In general, the trend is an increase 
in the spread rate with an increase in the tray width, presumably due to 
decrease in fractional heat loss to the side walls from the liquid and 
due to f^e viscous effects in the liquid at the side walls which alter 
the liquid phase energy convection process. When the flame temperature 
is low, e.g.: 30 °C, the pool width appears to be relatively incon- 

sequential except in the narrowest configurations. At higher fuel supply 
temperatures, however, the effect of tray width is felt over the entire 
range of our testing. Similar experiments by MacKinven [6] indicate that 
a tray has to be as wide as 60 cm before one can neglect the effects of 
the finite width on spread rate. 

Whereas the width of the tray produces marked change in the 
nature of the flames behind the spreading front, the leading edge of 
the flame appears to be early unaffected. (For narrow trays, the flame 
is laminar and its height is about 4-5 times the tray width. For the 10 
cm tray, the flame appears to be turbulent and its height is about 2-3 
times the tray width.) An estimation of the radiant heat transfer from 
the flames of different heights and of different characteristics is not 
made. 

Even though Figure 18 is prepared for a pool depth of 2.54 cm, 
the general features described above are valid for other depths as well. 

(h) Effect of Pool Depth : This is perhaps the most important 

effect which requires study both in theoretical and experimental view 
points. The reason for this importance is in the fact that in shallow 
pools the surface-tension-driven convection may become suppressed and 
hence the flame spread rate may be greatly reduced. Figure 19‘ shows the 
data cross plotted from Figures 15 to 17 at the narrowest width tested. 
Once again the points are not necessarily actual measurements but 
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but readings from the reference figures* 

When the pool is about 1 to 2 mm deep, ignition may be possible 
but flame spread is not possible except when the fuel temperature is 
near the flash point. MacKinven [61 points out the reasons for the ob- 
servation as (i) the excessive heat losses by conduction to the base of 
the tray and (ii) the suppression of the essential heat feed-forward 
mechanism of surface-tension-driven flow. 

At larger depths, the spread rates are correspondingly larger. 

For depths greater than about 6 mm, the flame spread rate increases with 
pool depth nearly linearly the slope being larger at higher fuel supply 
temperatures. Estimates show that the heat loss rate to the base of the 
pool becomes negligible for depths beyond about 1.5 cm. 

(i) Miscellaneous Observations : The flame spread experiment 

is conducted in trays of two different lengths, 122 cm and 61 cm. In 
almost all the cases, the fuel employed, namely butanol, permitted ignition 
by a simple exposure to a kitchen match flame. Based on the video records 
and visual observations, one can safely conclude that the ignition 
transient has extended no more than about 10-15 cm from the ignition end. 

In this length, there were observed liquid phase disturbances due to the 
preignition heat transfer. Once the flame is initiated and is spread 
over this length, a steady state is attained. All the flame spread rates 
reported above are taken from the flame location versus time data in 
the middle 1/3 of the length of the tray. As the flame would propagate 
over the length of the tray and approach' the far end of the pool, the 
end walls are expected to effect the surface-tensional convection and 
thus effect the flame spread. From an estimation of the dimension of the 
convection cell, we conclude that the last 10 cm or so of the pool length 
will result in a flame propagation altered by the presence of the end 
wall . 

6 • Observations in the Preignition Heating Experiments : As a 

propane flame is brought into the vicinity of a quiescent pool of liquid 
fuel, surface-tension-driven flow is initiated transiently. The ignition 
flame is supported by a controlled supply of propane through a capillary 
tube of i-d 3 mm. The resultant diffusion flame is about 2 cm long. The 
bottom part of this flame is nearly spherical in shape with a diameter of 
h out 5 mm. It is this nearly spherical region which acts as the energy 
source for the liquid fuel, the separation distance between the bottom of 
the flame and the liquid surface being about 5 mm also. Typical liquid 
phase flow patterns are already presented in Plate II. For all fuel 
supply temperatures tested (in the range of about 28 °C - 36 °C) the 
ignition time is in the range of 0.1 to 1 second* In fact, this fractionally 
variant ignition time is difficult to be measured accurately. Conse- 
quently, we satisfied outselves with visual and photographic observations 
of the development of surface-tension flows. 

Inasmuch as this compromise is made to limit our ignition 
experiments to observations of the liquid phase flow evolvement with 
time, we wondered if a heat source adjacent to a pool of water would pro- 
duce a transient surf ace-convection- induced flow similar to that induced 



in butanol. Our attempts failed to produce any distinct patterns but 
a similar unusual behavior of water has also been observed hy other 
investigators such as MacDonald [4], and Berg [45]. The inevitable 
presence of impurities in water which would greatly influence the sur- 
face tension is claimed by these authors to be responsible for the ob- 
served anamoly between theoretical predictions and experimental 
observations. 

7. Conclusions ; These experiments lead to the conclusions which are 
listed in the main body of the report (Chapter V). 
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